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ABSTRACT 


An improved version of the Wall Signature Method was developed to compute wall 
interference effects in three-dimensional subsonic wind tunnel testing of aircraft models in 
real-time. The method may be applied to a fullspan or a semispan model. 

A simplified singularity representation of the aircraft model is used. Fuselage, support 
system, propulsion simulator, and separation wake volume blockage effects axe represented 
by point sources and sinks. Lifting effects axe represented by semi-infinite line doublets. 
The singularity representation of the test article is combined with the measurement of wind 
tunnel test reference conditions, wall pressure, lift force, thrust force, pitching moment, 
rolling moment, and precomputed solutions of the subsonic potential equation to determine 
first order wall interference corrections. 

Second order wall interference corrections for pitching and rolling moment coefficient 
axe also determined. A new procedure is presented that estimates a rolling moment coef- 
ficient correction for wings with non-symmetric lift distribution. 

Experimental data obtained during the calibration of the Ames Bipod model support 
system and during tests of two semispan models mounted on an image plane in the NASA 
Ames 12ft Pressure Wind Tunnel axe used to demonstrate the application of the wall 
intexference correction method. 
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CHAPTER 1 


INTRODUCTION 

Wind tunnel tests have always played an important role in the development of mod- 
em aircraft. These tests are used to simulate atmospheric conditions experienced by an 
aircraft or spacecraft in free-flight. Aerodynamic forces and moments are measured and 
related to corresponding free flight values using Mach and Reynolds numbers. These mea- 
surements provide valuable information about expected performance, stability, and control 
characteristics of a new aircraft design. 

Large wind tunnel models, i.e. wing span on the order of 80% of the wind tunnel 
width, are often preferred in order to achieve a good simulation of viscous phenomena of 
the flow field. In this case, however, the presence of the wind tunnel wall and model support 
system change the free-air flow field experienced by the aircraft model. These flow field 
interference effects have to be considered to allow a reasonable comparison between wind 
tunnel test and free flight condition. Therefore, interference corrections to Mach number, 
dynamic pressure, and angle of attack have to be determined to improve test data quality. 

In the 1970s and 1980s, techniques were developed that use boundary measurements 
during a wind tunnel test to predict wall interference corrections. The Wall Signature 
Method introduced by Hackett et al. [1],[2],[3] and the Two-Variable Method introduced 
by Ashill [4], [5] were used extensively in 3-dimensional wind tunnel testing. Ashill [5] gives 
a detailed discussion and comparison of these techniques. 

The Wall Signature Method and Two-Variable Method are both based on potential 
flow theory. Computed wall interference corrections agree if each method is applied cor- 
rectly. However, a few differences exist between these two methods. Each method has its 
advantages and disadvantages. Table 1 compares important features of the Wall Signature 
Method and the Two-Variable Method. 

The Two-Variable Method does not require a singularity representation of the wind 
tunnel model to determine wall interference corrections. However, the wall interference 
correction calculation depends on an integration of the measured and interpolated surface 
pressure distribution on the wind tunnel wall. 
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The Wall Signature Method uses a singularity representation of test article, wall pres- 
sure measurements, precomputed perturbation velocities, and the principle of superpo- 
sition to compute wall interference correction. No integration and interpolation of the 
surface pressure distribution is required. Therefore, the Wall Signature Method is more 
flexible in the selection of wall pressure port measurements required for the wall interfer- 
ence calculation. Computed corrections are relatively insensitive to the number and the 
location of wall pressure ports used for the least squares fit. 

Table 1 : Comparison Wall Signature Method / Two— Variable Method 



Wall Signature 
Method 
( Hackett et al.) 

Two- Variable 
Method 
( Ashill et al.) 

Singularity Representation 
of the Wind Tunnel Model 

YES 

NO 

Measurement of the Flow Velocity 
at the Wind Tunnel Wall 

YES 

YES 

Number of Wall Pressure 
Port Measurements 

> 2 

> 100 

Wall Interference Correction 
Calculation based on a Perturbation 
Velocity Flow Field and the 
Principle of Superposition 

YES 

NO 

Wall Interference Correction 
Calculation based on an Integration 
of the Surface Pressure Distribution 
at the Wind Tunnel Wall 

NO 

YES 

Sensitivity of Computed Corrections 
to the Number and the Location of 
Wall Pressure Port Measurements 

LOW 

HIGH 


The Wall Signature Method was selected for use in the real-time Wall Interference 
Correction System (WICS) of the NASA Ames 12ft Pressure Wind Tunnel, because it is 
fast and does not depend on tin integration of a measured surface pressure distribution. 
The Wall Signature method will still provide corrections even if a large number of wall 
pressure ports do not provide a useful measurement. 

In general, the Wall Signature Method computes wind tunnel wall interference cor- 
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rections by introducing a simplified representation of the test article expressed in terms 
of singularities. Sources and sinks represent the fuselage volume and viscous separation 
wake blockage effects and horseshoe vortices or line doublets represent the lifting effects. 
In addition, power simulator blockage effects can be represented by a sink [6] . This sin- 
gularity representation is combined with a least squares fit of wall pressure measurements, 
data from calibration tests, and solutions of the subsonic potential equation, in the form of 
normalized perturbation velocities, to predict Mach number, dynamic pressure, and angle 
of attack corrections. 

During the past decade significant advances in the development of low-order panel 
method codes and computer hardware have made a fast calculation of complex three- 
dimensional internal flow field problems on workstation type computers possible. Panel 
method codes allow application of the Laplace Equation to realistic three-dimensional wind 
tunnel geometries which is important if the methodology of the Wall Signature Method is 
applied to the quasi-octogonal cross-section of the 12ft Pressure Wind Tunnel (PWT) at 
NASA Ames Research Center. It was shown by Ulbrich and Steinle [7], [8] that normalized 
panel method code solutions of the wind tunnel flow field combined with the Wall Signature 
Method can be used to predict subsonic wall interference corrections close to real-time. 

The revised formulation of the Wall Signature Method developed for the 12ft PWT is 
described in detail in this report. Figure 1 shows principle elements of the modified Wall 
Signature Method. Improvements of the Wall Signature Method were introduced to allow 
an application of the Wall Signature Method in real-time and to deal efficiently with a 
wide range of model and support system geometries. 

Originally, Hackett et al. [1] based their formulation of the Wall Signature Method on 
a “local” least squares fit procedure. They introduced a piecewise approximation of the 
wall signature using a parabola for its maximum and a tanh - function for its downstream 
asymptote. The location of singularities was found by matching the location of the maxi- 
mum of the parabola with the inflection point of the tanh - function. Unfortunately, this 
feature of the original formulation of the Wall Signature Method is difficult to use in a 
real-time correction system, as it requires the selection of wall pressure ports used for the 
“local” least squares fit of the maximum of the real-time wall signature. 

Ulbrich [9] introduced improvements to the Wall Signature Method to overcome the 
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limitations of a “local” least squares fit of the wall pressure signature. He suggested a 
“global” least squares fit procedure which matches the wall signature on all wall pressure 
ports using panel method code solutions of singularities placed inside the wind tunnel test 
section. In his approach, a “best” singularity location is found by minimizing the standard 
deviation of the least squares fit of the wall signature as a function of the singularity 
location. 

Support system wall interference corrections for fullspan model tests can also be found 
by applying the Wall Signature Method. In this case the Wall Signature Method has to be 
applied to the difference between the support system and the empty tunnel calibration at 
the wall pressure ports. Support system wall interference effects can be computed off-line 
and stored in a database. 

In the first part of this report, basic relationships of the proposed Wall Signature 
Method are derived for a fullspan and a semispan model. 

The second part of this report discusses the integration of the method into a wind tun- 
nel facility. Experimental data, obtained during tests of two different size semispan models 
mounted on an image plane in the NASA Ames 12ft Pressure Wind Tunnel (PWT), are 
applied to the modified Wall Signature Method to verify computed corrections. Exper- 
imental data recorded during the calibration of the Ames Bipod are also applied to the 
method. 
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CHAPTER 2 


WALL INTERFERENCE CORRECTION PREDICTION 
2.1 Definition of Interference Correction 


Wind tunnel tests allow the prediction of aerodynamic forces and moments acting 
on an aircraft model in atmospheric free-flight. Unfortunately, the wind tunnel wall and 
the model support system change the flow field experienced by the aircraft. Many of 
these changes can be ignored if the aircraft model is small compared to the wind tunnel 
height and width. However, if the span of the test article is large or if substanial flow 
separation occurs, wall and model support system interference effects cannot be neglected. 
Then, reliable estimates of interference corrections to Mach number, dynamic pressure, 
and angle of attack are necessary so that wind tunnel test data may be compared with 
free— flight conditions. 

In general, wall and support system interference corrections are defined as the dif- 
ference between the wind tunnel flow field and the free-air flow field experienced by the 
model (see Fig.2a) . Corrections are described in terms of a blockage factor c and an angle 
of attack correction a,- . Mach number and dynamic pressure corrections are related to the 
blockage factor computed at some reference point in the wind tunnel. For more detail on 
classical subsonic wall interference corrections, see AGARDograph 109, [10] . The block- 
age correction relates the free-stream velocity to a calibrated empty tunnel velocity 
U e at a model reference point v (see Fig. 2b) . The calibrated empty tunnel velocity U e 
captures the effects of the wind tunnel wall boundary layer growth, wall divergence, and 
orifice error. It is still necessary to correct for the wall interference effect of the test article, 
its separation wake, and the influence of the support system. 

The ratio between free-stream velocity Uoo and the calibrated empty tunnel velocity 
U e is expressed as a function of the blockage factor e, [10] : 


u e { v) 


1 + 


u «0) 

U e (u) 


1 + e(i/) 


( 1 ) 


where U{ is the axial velocity correction at the model reference point v caused by the model 
and support system interference effects relative to the calibrated empty tunnel velocity U e . 
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For small changes in velocity, second order approximations of Mach number and dy- 
namic pressure correction can be expressed as a function of the blockage factor e(i/) using 
a Taylor series expansion. These second order approximations are used if a large blockage 
factor is expected, e.g. during high angle of attack tests of aircraft models. We get (see 
Appendix 1) : 


Moo - Mjy) 
M e ( v) 


1 + 


7 - 1 


M\{v) 


<») + 


3(7 - l)Mg (i/) 


t{uf 


(2a) 


goo - ge(y) 

geM 


Ml{u) 


e(u) + 


1 - 


-M e V) + 


2-7 


Mt{u) 


( 2b ) 


The calculation of wall interference corrections based on the ideas of the Wall Signature 
Method requires the formulation of a blockage factor e and an angle of attack correction a t - 
such that a direct connection between experiment and the panel method code computation 
is possible [7], [8]. 

Studies by the author have shown that differences of the panel method code solu- 
tions can be used in combination with the Wall Signature Method to predict interference 
corrections. Figure 2b shows the relationship between the empty tunnel calibration, sup- 
port system calibration, wind tunnel test, free— air flow field, and the corresponding panel 
method code calculations using a simplified representation of the test article and support 
system in terms of singularities. This representation is uniquely defined if type, location, 
and strength of the singularities is known. The type and location of singularities must be 
specified by a test engineer. The Wall Signature Method is used to compute the strength 
of the singularities. The Principle of Superposition and panel method code solutions of 
the wall interference flow field axe used to determine wall interference corrections. 

It is assumed that singularities of the test article and support system are located inside 
a tunnel of constant cross-sectional area as effects of wall divergence and boundary layer 
growth are already included in the calibrated empty tunnel velocity U e . The calibrated 
velocity U e corresponds to a constant reference velocity U* of a constant cross-section 
wind tunnel. 

Comparing the flow fields depicted in Fig. 2b we can make the following approximation 
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of the velocity ratio U 00 (u)/U e (i') defined in Eq. (1) : 

„ u: + [ ( u: + x t(v ) ) - ( u; + O) ) ] 


UM 


u: 


= i + 


<(") - «m0) 


(3a) 


u: 


The total perturbation velocity u* of the wind tunnel flow field caused by the test article 
and support system is expressed as the sum of the perturbation velocity contribution u* m 
of the test article and the perturbation velocity contribution u* s of the support system (see 
Fig. 2c) . The perturbation velocity component uj, of the support system perturbation 
velocity flow field can further be represented as the sum of the perturbation velocity 
component — u* due to the wall interference of the support system and the perturbation 
velocity component u* due to the direct influence of the support system (see Fig. 2d). 
Finally we can write : 

O) - o) + + «:(*) i 


u oo(u) 
UeW 


= 1 + 

= 1 + 


u: 


r - U* m (v) 


u; 


+ 


<(») - 


u: 


+ 


u: 


(36) 


It is difficult to predict the direct influence of the support system on the model flow 
field in terms of a perturbation velocity component u* . The Wall Signature Method 
cannot be used to determine the direct influence of the support system on the model as 
this technique has been developed for wind tunnel wall interference studies. Experimental 
or CFD studies have to be used to estimate the direct influence u* . 

For the present study it is assumed that the direct interference between support system 
and wind tunnel model is small, i.e. u* ft; 0 . Then, comparing Eqs. (1) and (3b), we get 
for e(v) : 


e(u) 


r <m0) - O) i 

\ 

U* s (i/) - <( 1/) j 

l u: \ 


u; J 


= e m (j/) + e $ (u) (4) 


Interference velocity components in pitch and yaw axis direction can be expressed in 
forms similar to Eq. (4) if direct influence of the support system is neglected, i.e., v* ft; 0 
and w* ft: 0 . We get : 


Ue(u) 


- KM 
u: 


M 


u: 


(5a) 
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w. 


M 


Ve( v ) 


KM - KM i , r KM - w M 


u: 


+ 


u: 


= a m + a. 


(56) 


The velocity ratios a m and a { are angle of attack corrections due to the model and support 
system wall interference effects. 

The interference velocity component in the z-axis direction is related to the free- 
stream angle of attack aoo experienced by the test article at the model reference point. 
Therefore we get : 

otoo(v) = a t ( i/) + (6) 


where a t is the geometric angle of attack measured relative to the wind tunnel centerline 
and dj is the angle of attack correction due to lift interference of the test article and due to 
the change of the flow field angle caused by the interference flow field of the model support 
system. The total angle of attack correction of a fullspan model is then : 


a M 


Wi(u) 

Uoo(v) 


U e (i/) Wi(u) _ U e (v) r 

Uooiv) ' UM = u„( v ) ' 1 + a * 


( 7 ) 


where U e (i/)/Uoo(v) and Wi(i/)/U e (u) are given by Eqs. (3b), (5b) . The angle of attack 
correction for the semispan model is obtained by replacing Wi(v) by u t -(i/) in Eq. (7) . 

The calculation of the Mach number, dynamic pressure, and angle of attack correction 
using Eqs. (2a), (2b), (4), (7) is reduced to finding blockage corrections e m and e, and angle 
of attack corrections a m and a s . 

Model corrections e m and a m , i.e. velocity differences [uj m — u^/U*, [v t * m — v^l/U* 
or [u/f m — w^l/U*, and support system corrections e, and a 4 , i.e. velocity differences 
— u*]/U*, — v’]/U* or [to*, — w*\/U* , can be computed by using the modified Wall 

Signature Method. This is possible because the wall interference flow field caused by the 
wind tunnel model and support system can be treated as a far field effect. This will be 
explained in detail in the following sections. 


2.2 Panel Method Code Solution 


In the previous section, it was demonstrated that wall interference effects of a test arti- 
cle and support system , i.e. e m , e s , a m , and a,, can be found by calculating dimensionless 
velocities [u* m - K m - <,]/[/;, [w* tm - w^]/U;, [u t % - <]/f/ e *, [t> t * f - v*]/U*, 
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and [wj, — w*]/U* at a selected test article reference point v. These dimensionless veloc- 
ities axe computed by superimposing panel method code solutions and applying the Wall 
Signature Method. In real-time operation, the Wall Signature Method uses a singularity 
representation of the test article in combination with the measurement of wall pressure, lift 
force, propulsion simulator thrust force, pitching moment, and precalculated normalized 
perturbation velocities to predict model wall interference corrections. The Wall Signature 
Method may also be used to predict support system wall interference corrections by talcing 
the difference between the wall pressure port calibration of the support system and empty 
tunnel. However, the application of the Wall Signature Method is only possible, if precal- 
culated solutions of the subsonic potential equation in the form of normalized perturbation 
velocities are linear with respect to singularity strength. 

In general, the perturbation flow field of a singularity placed inside a wind tunnel 
of constant cross-sectional area is a linear function of the singularity strength. Figure 3 
depicts a singularity of strength cr 0 located inside a selected wind tunnel configuration of 
constant cross-sectional area A . The corresponding flow field solution is obtained by using 
a modified panel method code as a boundary value problem solver (see Appendix 2 for a 
detailed description of these modifications). Figure 3 also shows a singularity of the same 
type and location but with different singularity strength oq . The strength a 0 and <j\ of 
these two singularities are related as follows : 


( 7 1 = A • <t 0 (8a) 

We know for the corresponding panel method code solutions in terms of the perturbation 
velocities : 

u* = \-u * 0 (86) 

Equation (8b) is valid as long as the cross-section of the selected panel geometry of the wind 
tunnel test section is constant. This relationship allows the user to find all wind tunnel 
flow field solutions of <J\ 7^ 1.0 [ft 2 /sec] by simply multiplying the flow field solution of 
(T 0 = 1.0 [ft 3 /sec] with the factor A = cr\/cr 0 . 

In the following section, details will be developed as to how the Wall Signature 
Method can be combined with the linear property of panel method solutions defined by 
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Eqs. (8a), (8b) to obtain test article and support system wall interference corrections for a 
a fullspan or semispan model. 


2.3 Fullspan Model 

2.3.1 Test Article Wall Interference Correction 

The Wall Signature Method can be applied to compute the blockage factor e m and 
angle of attack correction a m caused by test article wall interference effects at a reference 
point V* . 

The Wall Signature Method uses a simplified representation of the test article in 
terms of point sources, point sinks, and line doublets in combination with measurement of 
wall pressure, lift force, propulsion simulator thrust, pitching moment, and precalculated 
normalized perturbation velocities to obtain the blockage factor and angle of attack cor- 
rection. The blockage factor, e m , and angle of attack correction, a m , are linear functions 
of the singularity representation of the test article (see previous section). The principle of 
superposition allows e m and a m to be expressed as the sum of contributions of “n” singu- 
larities that represent the test article. Figure 4 shows, as an example, how the principle 
of superposition can be applied to a simplified representation of a test article using three 
singularities. Knowing that e m and a m are a linear function of the singularity strength <r, 
we get for singularities : 

n 

f m(^) = ^2 (Tk • W(l', fc) (9a) 

jt=l 

Tl 

<*m(^) = (Tk ■ fc) (96) 

Jfc=l 

where ul(v, k) and wi(v, k ) Eire normalized perturbation velocities of the wall interference 
flow field. These normalized velocities are dimensionless; perturbation velocities divided by 
unit singularity strength per unit velocity. Perturbation velocity ul(v, k ) is the normalized 
axial perturbation velocity component of the wall interference flow field and k) is the 
normalized perturbation velocity component perpendicular to the wing plane of the test 
article. Normalized perturbation velocities have the unit [1/m 2 ] or [l//< 2 ] . Singularity 
strength <rjt has the unit [m 2 ] or [ft 2 ] and is compatible with these normalized perturbation 
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velocities. Singularity strength cr* is defined as singularity strength cr in [m 3 /sec] or 
[ft 3 / sec] (see also App. 3,4) divided by reference velocity U* or U re j '• 

cr cr 

° k = U* = U^l 

The strength of the singularities is computed similar to the procedures outlined by 
Ulbrich and Steinle [7], [8], and Koning [6], using the measurement of lift force, thrust 
force, pitching moment, and a least squares fit of wall pressure signature. Figure 5a 
summarizes the basic steps in the application of the Wall Signature Method to a fullspan 
model configuration. Rectangular boxes in Fig. 5a symbolize real-time measurements, 
rounded boxes symbolize information stored in database files, and elliptical boxes symbolize 
computational procedures. 

Fuselage volume blockage effects are represented by point sources and point sinks. The 
location of these sources and sinks has been selected by the user such that pairs of sources 
and sinks axe related to Rankine bodies describing the fuselage volume of the test article. 
Strengths <r\, ...,cr^ of point sources and strengths cr v + 1 , ..., of point sinks axe reduced 
to a single variable cr* if weighting factors uq, ..., u>2*> are introduced. These weighting 
factors must be defined by the user. We then obtain : 

— J — — wj ; 1 < j < 2 tj ( 11 °) 

cr* 

where 

Wj = - Wj- n 5 Tj + l < j < 2rj (life) 

Separation wake blockage effects are represented by a set of sources, o^rj+i, ..., . 

Again, weighting factors axe introduced to reduce the number of independent variables. 
Assuming that the strength cr*, is a common reference strength of sources related to the 
separation wake we get : 



= wj ; 2r] + 1 <j < £ (12) 

cr** 

The weighting factors u^+i, ...,W{ axe, by definition, greater than zero since separation 
wake blockage effects are modeled as sources. The calculation of the strength of a total 
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number of point sources and sinks representing fuselage volume and wake blockage 
effects is then reduced to finding the values of cr* and cr** using the Wall Signature Method. 

A least squares fit of velocities derived from wall pressure measurements is used to 
calculate a * and <r„, (see Fig. 5b) . It is assumed that a total number of “n” singularities 
has been selected to represent the test article. Axial velocities are derived from pres- 
sure measurements at V wall pressure orifices. The measured real-time wall signature, 
U e (6) + itt(<5) , is corrected for orifice error, wall divergence, and wall boundary layer 
growth by subtracting the measured velocity, U e (8) -f u tJ (<5) , of the support system 
calibration at each wall orifice location “6”. The computational equivalent, U*, to the 
calibrated velocity, U e or U re f , at the wall pressure orifice “5” is constant everywhere 
inside of the test section as it is computed using a constant cross-section wind tunnel. 
Normalized perturbation velocities of the wind tunnel flow field at wall pressure orifice 
locations, “6” , have to be introduced to compute the strength of singularities representing 
the test article. They relate the strength of point sources, sinks, and fine doublets to the 
perturbation velocity components at the wall pressure orifices. The strengths, cr^+i , ..., cr^, 
of line doublets representing lifting effects of the test article are estimated by combining 
lift force, pitching moment, and rolling moment measurements with the Kutta/ Joukowski 
formula (see Appendix 5 and Appendix 19) . The strengths, 1 , ..., u n , of point sinks 
representing blockage effects of propulsion simulators axe estimated using thrust measure- 
ments and Koning’s formula (see Appendix 10) . The normal equation of the modified 
Unear least squares problem depicted in Fig. 5b is giver, as, [11] : 



(13a) 

(136) 

(13c) 

(13d) 

(13c) 
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( 13 /) 


b, 


Bmx 1 — 



[ U e (S) + U t (6) ] - [ U e (S) + u u (6) ] 

Uref 


y o-k-u t (6,k) 

k=Z + 1 


( 13 *) 


The vector X contains the strength of singularities modeling the fuselage volume and 
wake blockage effects. The matrix A contains normalized perturbation velocities of the 
wind tunnel flow field. The vector B contains residual perturbation velocity components 
caused by the fuselage volume and separation wake. Line doublet and propulsion simulator 
contributions, k = £ + 1, ..., n, are subtracted from the measured wall signature difference 
between the wind tunnel flow field and the support system calibration at the wall pressure 
ports “6” (see Fig. 5b). 

Perturbation velocity, ui{8, fc), is the normalized perturbation velocity of the wind 
tunnel flow field of singularity “jfe” at wail pressure orifice “<5” . It is defined as the 
dimensionless perturbation velocity divided by unit singularity strength per unit velocity 
at wail pressure orifice “<5” due to a singularity “it” located inside the test section. 

In general, it is required that the measured velocity, U e (6 ) + u t (6), at wail pressure 
port “5” can be approximated by its component in the streamwise direction. In practical 
applications, however, a least squares fit can tolerate a few wall pressure measurements 
that do not fulfill this condition as the normal equation of the least squares fit, Eq. (13a), 
assigns equal weight to all wall pressure port measurements. 

The solution of the two variable linear least squares problem defined in Eq. (13a) can 
be written in explicit form as : 


L 2x 1 


= 




[A T ' B] 


2x1 


(14) 


The solution vector X is computed by using the Singular Value Decomposition tech- 
nique [12], which is the numerical method of choice for linear least squares problems. 

It is now possible to compute the blockage correction, e m (is), and the angle of attack 
correction, a m (i/), at reference point “i/” as the strength and location of all singularities 
representing the test article and wake are known. 
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2.3.2 Support System Wall Interference Correction 


The Wall Signature Method can also be used to find the blockage factor, e„, and the 
angle of attack correction, a,, caused by support system wall interference effects. The 
application of the Wall Signature Method to the support system wall interference problem 
closely follows procedures discussed in the previous section. Figure 6a summarizes basic 
elements in the application of the Wall Signature Method to the support system wall 
interference problem. 

Support system wall interference effects have been defined in Eqs. (4) and (5b) . 
Similar to Eqs. (9a), (9b), e, and a, are a linear function of the singularity strength. It 
is assumed that a total number of sources and sinks of unknown strength axe used to 
represent blockage effects of the support system. We then get : 




«(* / ) 

= } (Tk • k ) 

k= 1 

(15a) 


= ^2 ak ■ *) 
Jfc=l 

(156) 


where u,( v, k ) and Wi(u, k ) are normalized perturbation velocities of the wall interference 
flow field. 

Again, positions of sources and sinks modeling solid volume blockage effects are se- 
lected such that pairs of sources and sinks are related to Rankine bodies describing the 
volume of the support system (see Eq. (11a), (lib)). Sources related to wake blockage 
effects of the support system axe placed where a wake separation on the support system is 
expected (see Eq. (12)). Introducing weighting factors it is possible to reduce the number 
of unknown singularity strength values to two, i.e. a* and <r*, . 

Assuming that “m” wall pressure measurements are taken during the support system 
calibration and that singularities are used to represent the support system, we get the 
following normal equation of the least squares problem, [11] : 

[■A. 2 xm ' A mX 2] ’X2xl = A^ > 2xm ' B mx i (16a) 


where 


B 


m x 1 



(166) 
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bs 


(16c) 


[ Ue(6) + u t ,(S) 1 - U'(6) 

U ref 

The vector X contains the strength of singularities describing the support system volume 
and wake blockage effects (see also Eq. 13b) . The matrix A is given by Eqs. (13c), (13d), 
and (13e) using the perturbation velocities ut(6, k) of the singularities 1 < k < £ represent- 
ing the support system. The vector B contains perturbation velocity components caused 
by support system volume and separation wake. These perturbation velocity components 
are the measured wall signature difference between the support system calibration and the 
empty tunnel calibration at wall pressure ports “<5” (see Fig. 6b). 

2.4 Semispan Model 

In general, blockage and angle of attack corrections are computed using Eqs. (4) and 
(5b) if the modified Wall Signature Method is applied to a semispan model. Studies by 
the author have shown that a semispan model mounted on a finite length image plane 
may be treated similar to the fullspan model configuration. It is only necessary to select 
the proper geometry of the wind tunnel, i.e. the cross-section of the wind tunnel channel 
above the image plane surface plus its reflected image, for the calculation of normalized 
perturbation velocities (see Fig. 7) . However, a new calibration of the empty tunnel 
velocities, U e (S) and U re j, must be conducted because the installation of the image plane 
changes the empty tunnel geometry. 

Semispan models normally consist of half of the fuselage mounted on the image plane. 
It is therefore necessary to place singularities representing the fuselage volume on the 
surface of the image plane. This requires further modification of a panel method code to 
compute normalized perturbation velocities (for more detail see Appendix 2) . 

No support system is present in the test section, i.e. u t = u tm , u ts = 0.0, and u, = 0.0 
(see boundary value problems depicted in Fig. 8) . Thus Eq. (4), (5b), (7) are replaced by 
the following expressions : 


e(v) ss 

[ - O) 1 

i v:( v) \ 

= fm(t') 

(17) 

Vi{v) 

V ' C-w- 

\ - O) i 

1 u;(v) J 

= a m (y) 

(18a) 
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ai(u) = 


Vi(u) 


u e (v) I ?i(t/) 


VM 


a m (v) (186) 


«oo(^) Uoo(v) U t {v) Uoo(i') 

Normalized perturbation velocities, tZ7, v7, and u7, of the interference and wind tunnel 
flow field are computed using a panel method code. These perturbation velocities axe 
required for the least squares fit of the wall signature and the calculation of wall interference 
corrections. The user has to make sure that the orientation angle of the line doublets is 
changed for semispan tests (see Appendix 4). 

Figure 9a summarizes the basic elements of the Wall Signature Method for semispan 
model tests. The least squares fit of the wall signature, defined in Eqs. (13a) to (13f), still 
applies. Again, it is necessary to subtract the contribution of the wing line doublets and 
the contribution of the propulsion simulator point sinks from the wall signature difference. 
However, the image plane is calibrated as a part of the wind tunnel wall and no support 
system is present in the tunnel. Therefore, the wall signature difference is defined as 
the difference between the wind tunnel measurement, U e (S ) + «t(6), and the image plane 
calibration, U e (6), (see Fig. 9b). Equation (13g) is replaced by 


bs 


[ U e (S) + u t (S ) ] - UejS) 

U ref 


T, &k ' U t (6, k) 

t=(+l 


(19) 


For a semispan model configuration the lift force measurement and the Kutta/Joukowski 
formula are needed to determine the strength of the line doublets of the wing (see Ap- 
pendix 5). 


2.5 Compressibility Effects 

Compressibility effects have to be taken into account if the modified formulation of 
the Wall Signature Method is applied to a subsonic wind tunnel flow field. 

In general, the Prandtl/Glauert transformation may be used to approximate the effects 
of compressibility in a flow field. This transformation requires a coordinate stretching of 
the wind tunnel geometry, singularity location, singularity strength, and reference point 
location when normalized perturbation velocities are being computed as outlined in the 
previous sections (for more detail on the transformation see Ref. [13] and [14]) . 

The application of the Prandtl/Glauert rule to the flow field of a singularity located 
inside of a wind tunnel requires several steps. First, coordinates (x,y,z) of the tunnel 


16 



geometry and of the singularity position are transformed to corresponding incompressible 
coordinates (x,y,z) using the Prandtl/Glauert transformation : 

x — x (20a) 

y = y • y/l - M 2 = y- (3 (206) 

z — z • \/l — M 2 = z • (20c) 

It is necessary to stretch the strength of the singularity as well. The strength of a 
source is related to the cross-sectional area of a corresponding halfbody, [13], as 

a = IT • R 2 ■ Uoo = K ■ [ y 2 + Z 2 ] • Uoo (21) 

where R equals the radius of a halfbody far downstream of the source location. Singularity 
strength a has the units [m 3 / sec] or [ft 3 /sec] . Comparing Eqs. (20b), (20c), (21) we see 
that the singularity strength is proportional to ft 2 . The strength of a line doublet is 
related to a corresponding wing span increment, i.e. Ay, and angle of attack, i.e. w/Uoo 
or d z/d x . Considering Eqs. (20b) and (20c) we see that the line doublet strength has to 
be proportional to (3 2 . Finally, we get for the incompressible singularity strength : 

d = a • [ 1 - M 2 ] = <j ■ j3 2 (22) 

A Mach number, M, derived from the calibrated velocity, U re j, at a test section 
reference point should be used for the coordinate and singularity strength stretching. 

Next, the incompressible flow field is obtained by using the modified version of the 
panel method code PMARC,[15] as a boundary value problem solver. Finally, perturbation 
velocities calculated relative to the inlet velocity are transformed from the incompressible 
flow field, (u,w,it)), to the compressible flow field, («,«,«;) . We then get : 

(23a) 

(236) 
(23c) 


P 2 


v 



w 
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A quantitative check of compressibility effects is possible if we place a point source 
of unit strength inside the geometry of the 12ft PWT and use a modified panel method 
code and the Prandtl/Glauert transformation described above to solve the flow field. Fax 
downstream of the point source the flow field is essentially one-dimensional and asymptotes 
of perturbation velocities for different Mach numbers can be calculated. The results of the 
panel method code calculations on Row 1 (see Fig. 17b for location of Row 1) for Mach 
numbers M = 0.0, 0.3, 0.6 are plotted in Fig. 10 . Asymptotic values of perturbation 
velocities can be found in Table 2 . 

Independent of the panel method code calculation it is also possible to apply the 
Area- Velocity-Relation of compressible flow, [16] . We can write : 

d u 1 — d A , 

Uoo = 1 - M 2 ' A 

The 12ft PWT has a cross-sectional area A = 109.74 [ft 2 ] . A point source of strength 
a /Uoo = 1.0 [ft 2 ] corresponds to a change of cross-sectional area of —dA = 1.0 [ft 2 ] (see 
also Eq. (21)). Equation (24) is applied to Mach numbers M = 0.0, 0.3, 0.6 . Results axe 
shown in the table below. 


Table 2 : Calculated Asymptotes of Perturbation Velocity 


Mach Number 

Panel Method Code 

Equation (24) ; du/Uoo 

0.0 


0.0091 





0.0140 

0.0142 


Comparing the results of the panel method code calculation and the Area- Velocity- 
Relation, we get excellent agreement to verify application of the Prandtl/Glauert trans- 
formation. 

2.6 Normalized Perturbation Velocity Definition 

In general, normalized perturbation velocities can be defined as dimensionless pertur- 
bation velocities divided by unit singularity strength per unit velocity caused by some sin- 
gularity. For example, the normalized perturbation velocity, uj(6, i fc), at a point, (x«, ys , zs), 
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in the free-air flow field of a unit strength source located at (xk,yk, z k) can be defined 
as (see also Eq. (3.2a) in Appendix 3) : 


u s (6,k) 


u(S, k) Uref _ ujSJe) 1 

UreJ O' Uref O k 

1 X$ — Xk 


4 * [[x s -xk] 2 + [ ys-yk ] 2 + [z« - zt] 2 ] 3/2 


(25) 


Two different sets of normalized perturbation velocities for a given singularity type, 
location, and specified model reference point location have to be computed if the Wall 
Signature Method is applied. The first set is related to the wind tunnel flow field of a 
singularity. The second set is related to the wall interference flow field of a singularity. 

Normalized perturbation velocities, Ut(6, k), of the wind tunnel flow field of a singu- 
larity axe required for the least squares fitting of the wall signature on wall pressure orifices 
“ 6 ” (see Eqs. (13d), (13e), (13g), (19)) . Each normalized perturbation velocity can be 
interpreted as a dimensionless streamwise perturbation velocity divided by unit singularity 
strength per unit velocity at wall pressure orifice “ S n caused by the wind tunnel flow field 
of a singularity located at position “ife” . 

Normalized perturbation velocities, Ui(v, k ) and wl(v, k), of the wall interference flow 
field are required to determine wall interference corrections at reference point “i/” (see 
Eqs. (9a), (9b), (15a), (15b)) . Each of these perturbation velocities can be interpreted 
as a dimensionless perturbation velocity component divided by unit singularity strength 
per unit velocity at flow field reference point “i/” caused by the wall interference flow field 
of a singularity located at position “ife” . Normalized perturbation velocity ui(v, k ) is the 
streamwise perturbation velocity component and normalized perturbation velocity Wi(v, k ) 
is the perturbation velocity component perpendicular to the wing plane. 

Explicit equations of these normalized perturbation velocities can be derived in the 
case of a wind tunnel with rectangular cross-section using the Method of Images. However, 
normalized perturbation velocities of a wind tunnel with non-rect angular cross-section 
have to be computed using a panel method code as a boundary value problem solver (see 
Ref. [15], [17], and Appendix 2 for more detail) . 

Exact Bessel Function solutions of the angle of attack correction of a line doublet 
located in plane z = 0.0 inside of a wind tunnel with circular cross-section are available. 
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NACA TN 2454, [18] lists these corrections in the form of upwash factor tables. Tables 
can be compared with the normalized perturbation velocities of the wall interference flow 
field that were computed using a modified panel method code [15], [17] . 

At first, an upwash factor table used for fullspan model tests (Table I on p.32 of NACA 
TN 2454, [18]) is compared with the corresponding panel method code solution. Figs. 11 a 
and 1 lb compare upwash factor F as a function of the dimensionless streamwise coordinate 
and the dimensionless line doublet location (lateral coordinate 77 = 0.7). The exact Bessel 
Function solution and the numerical panel method code solution of the upwash factor F 
show reasonable agreement verifying the normalized perturbation velocity definition of the 
fullspan configuration. 

Table 3 compares the input and accuracy characteristics of NACA TN 2454 and the 
panel method code solution of the normalized perturbation velocities. 

Table 3 : Comparison NACA TN 2454 / Panel Method Code 



NACA TN 2454 

Panel Method Code 

Tunnel 

Geometry 

Circular Tunnel 
One Bipolar Tunnel 

Any Tunnel Geometry 

Solution 

Exact (Circular) 

Numerical 

Type 

Approximation (Bipolar) 

Solution 

Singularity 

Type 

Line Doublet 

Point Source, Point Doublet, 
Line Doublet 

Singularity 

Z=0.0 (Circular) 

Minimum distance from wall 

Location 

Y=0.0 (Bipolar) 

panels ss 0.2 x tunnel radius 

Reference Point 

Z= 0.0 (Circular) 

Minimum distance from wall 

Location 

Y= 0.0 (Bipolar) 

panels ~ 0.2 x tunnel radius 


Unfortunately, no rigorous solution of the upwash factor F, i.e. angle of attack cor- 
rection, of a line doublet located inside of a bipolar wind tunnel (image plane / semispan 
configuration) is available. However, NACA TN 2454 provides an approximation of the 
upwash factor F for a bipolar wind tunnel that can be compared with the results obtained 
by applying a panel method code. Fig. 11 c compares this approximation of upwash factor 
F (Fig. 5(c) on p.56 of NACA TN 2454, [18]) with the panel method code solution. Both 
approximations show reasonable agreement verifying the normalized perturbation velocity 
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definition of the semispan model configuration. 

A panel method code allows the user to compute normalized perturbation velocities for 
any type of constant cross-section wind tunnel geometry that can be paneled. Singularities 
and reference points can also be placed anywhere inside of the wind tunnel as long as the 
minimum distance of the singularity or reference point from the paneled wind tunnel wall 
(60 panels used to represent tunnel cross-section) is greater than 0.2 x tunnel radius. 
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CHAPTER 3 


APPLICATION OF THE METHOD TO WIND TUNNEL TESTS 
3.1 Real-Time Wall Interference Calculation 

The revised and improved version of the Wall Signature Method presented in this 
report can be used to predict the Mach number, dynamic pressure, and angle of attack 
correction at a test article reference point due to the subsonic wind tunnel wall interference 
effects. Post-test analysis of subsonic wall interference effects is also possible as long as 
the exact position of the test article in the wind tunnel is known for a specific angle of 
attack setting. The location of singularities representing the test article is directly related 
to this position. Perturbation velocities of the wind tunnel and interference flow field can 
be computed as outlined in the previous chapters. 

The real-time calculation of wall interference corrections is fast because the Wall 
Signature Method only requires superposition of the perturbation velocities, application 
of the Kuita/Joukowski formula (Appendix 5), Koning’s formula (Appendix 10), and the 
solution of a 2 x 2 linear system of equations related to the least squares fit of the wall 
signature. 

The precalculation of normalized perturbation velocities used for the real-time least 
squares fit of the wall signatures and for the calculation of corrections has to be done on a 
mainframe computer or fast workstation since a realistic implementation of the proposed 
Wall Signature Method requires the calculation of perturbation velocities for many different 
singularity types, locations, and Mach numbers. Figures 12a, 12b depict geometries of the 
NASA 12ft Pressure Wind Tunnel test section that are selected for the calculation of the 
perturbation velocities of the wind tunnel and interference flow field using a panel method 
code as a boundary value problem solver. Precomputed normalized perturbation velocities 
have to be stored in a database that is accessed during a wind tunnel test. 

A singularity and reference point grid has to be used for the calculation of the per- 
turbation velocity database. These two grids should be selected such that they allow for 
a real-time interpolation of all conceivable singularity and reference point locations in the 
wind tunnel test section. Perturbation velocities required for singularities representing the 
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test article change as a function of test article geometry and position in the wind tun- 
nel test section. Size, complexity, and accuracy requirements of the perturbation velocity 
database have to be balanced to guarantee best real-time performance. Thus, perturbation 
velocities of the wind tunnel and interference flow field for a given test article location are 
found in real-time by applying a tri-linear interpolation (real-time singularity position; 
Appendix 6) and parabolic interpolation (real-time Mach number; Appendix 7) to the 
precomputed perturbation velocity database. 

Blockage effects of the support system can be computed off-line by applying the Wall 
Signature Method to the difference between the support system and the empty t unn el 
calibration. Computed support system wall interference corrections on the reference point 
grid have to be stored in a database as a function of the support system calibration 
variables. The support system wall interference corrections are added in reed-time to the 
wall interference corrections caused by the test article. 

Post-test analysis of the interference effects is based on minimizing the standard 
deviation of the least squares fit of the wall signature as a function of the location of the 
test article singularities (see Appendix 8). This procedure provides an optimal singularity 
representation of the test article. 

Studies of the author have shown that the real-time speed of the wall interference 
calculation is governed by the efficiency of the interpolation of the perturbation velocities 
using the precomputed perturbation velocity database. 

Figures 13a, 13b depict basic elements of the real-time Wall Signature Method for 
a fullspan and a semispan model configuration. The empty tunnel calibration, support 
system calibration, read— time wall pressures , lift force, propulsion simulator thrust, and 
pitching moment measurements are the critical link between experiment and the panel 
method code calculation. Matching conditions between the wind tunnel test and the panel 
method code calculation will be discussed in detail in the following section. 

3.2 Matching Conditions 

The least squares fit proposed for a fullspan and semispam model configuration relates 
experimental data, i.e. wadi pressure measurements, to precomputed normalized perturba- 
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tion velocities. Therefore suitable equations have to be found to convert the wall pressure 
measurements to perturbation velocities. 

Perturbation velocity differences [U e + u t ] — [U R + ut $ \ > \U e + uu\ — U e , and [t/ e + u ( ] — 
U e defined in Eqs. (13g), (16c), and (19) can be related to the wall pressure measurements 
taken during calibration and real-time wind tunnel test by applying the energy equation 
and the isentropic flow assumption. Assuming that the total temperature Tt and total 
pressure pr in the tunnel settling chamber, and the static pressure p at wall pressure orifice 
6 axe known, we get for the flow velocity (see Fig. 14a) : 


U(p(6)) 


2 j R Tt 
7-1 


1 - 


1 - 


PT - p(5) 


PT 


(7-l)/7 


(26) 


Dimensionless perturbation velocities can then be written as : 


[ U e (6 ) + «t(5) ] - [ U e (S) 

+ ] 

U(ptun{6)) - U(p sup (6)) 

(27a) 

U ref 
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[ U e (6 ) + ut s (6 ) ] - 

■ U e (6) 

U(psup($)) — U(p em p(6)) 

(276) 

U ref 


U ref 

[ U e (6) + u t {6) } - 

Ue(6) 

U( Ptu n (6)) - U( Pernp (6 )) 

(27c) 

U ref 


U ref 


The flow velocity U(p(6)) in Eq. (26) is written as a function of the pressure difference 
pr — p(S) as it is easier to measure a pressure difference at a wall pressure port. Real-time 
static pressure p ttm , support system static pressure p sup , and empty tunnel calibration 
static pressure p em p are recorded at each wall pressure port “6” . Measured velocity 
U e {S) = U (pemp (<5)) at a wall pressure port “5” is the first matching condition. It is 
required to obtain the perturbation velocities defined in Eqs. (27b), (27c) . 

The velocity U re j in Eqs. (13g), (16c), and (19) is the second matching condition 
between the wind tunnel flow field and the corresponding dimensionless panel method 
code calculation. It is required to non-dimensionalize the perturbation velocities defined 
in Eqs. (27a), (27b), (27c) . It can be considered as the constant flow velocity inside of 
a hypothetical constant cross-sectional test section. This velocity should be measured 
during the empty tunnel calibration at a specific model reference point 9 (see Fig. 14b; 
this point has the coordinates X=120.71 [ft], Y=0.0 [ft], Z=0.0 [ft] in the 12ft PWT). The 
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velocity U re f can be measured by using, e.g., a static pipe installed in the wind tunnel 
during the calibration. We then get : 


Uref 


2 7 R Tp 
7-1 


1 


1 


PT Pernp(9) 


(7-l)/7 


PT 


(28) 


The real-time lift force and pitching moment measurements are related to the strength 
of the corresponding line doublets using the Kuita/Jouhowski formula. Its application re- 
quires knowledge of velocity U re f and density p re f (9) at a model reference point 9 measured 
during the calibration of the wind tunnel. This velocity and fluid density connect lift force 
and pitching moment measured in [lbf],[ft * Ibf] or [A r ],[A r * m] to the definition of the 
line doublet strength in [ft 2 ] or [m 2 ] (see Appendix 5). The fluid density is also required 
to relate the propulsion simulator thrust measurement to the sink strength if Koning’s 
formula is applied (see Appendix 10). The density p rc /( 9) is the third matching condition. 
It is found by applying the ideal gas and the isentropic flow relationship at a test section 
reference station 9 . We get : 


Pref (9) — 


PT 

R T t 


Pemp(9) 


1/7 


PT 


(29) 


where p em p(9) is the static pressure measured at the model reference point 9 . 

It is interesting to note that the perturbation velocity differences [t/ e (<5) + ut(<5)] — 
[C/ e (5) + ittj(^)] 5 [f/ e ( 6 ) + Uts(£)] — U e (6) , and [I/ e (5) + ut(6)] — U e (6) remove the influence 
of the wall boundary layer growth, orifice error, image plane, and wadi divergence from the 
least squares fit. Therefore it is possible to use the geometry of an equivalent wind tunnel 
with constant cross-sectional area for the calculation of normalized perturbation velocities 
of the wind tunnel and interference flow field. 

Three matching conditions, i.e. U e {6), U re f, and p re f , establish a link between the 
measurement of wall pressures, forces, moments and the panel code solutions of the wind 
tunnel and interference flow field expressed as normalized perturbation velocities. Figure 
14c summarizes the importance of these matching conditions. 


3.3 Application of the Method to Semispan Models 


In the summer of 1996, two different sized semispan models were tested in the NASA 
Ames 12ft Pressure Wind Tunnel (PWT). Both models, i.e. the 8 % and 14 % scale 7J7 
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semispan models, were provided by the Boeing Corporation. Each model was mounted on 
the image plane in the 12ft PWT. Figure 15 shows a similar test configuration. 

Both models were tested over a wide range of angle of attack, total pressure, and Mach 
number settings. For the present study, two runs were selected. During Run No. 154 the 
8 % scale model was tested from —20.28° to 19.82° uncorrected angle of attack at a total 
pressure of 2.0 [atm] and a Mach number of 0.25 . During Run No. 219 the 14 % scale 
model was tested from —4.03° to 9.98° uncorrected angle of attack at a total pressure of 
2.0 [atm] and a Mach number of 0.30 . 

The application of the present wall interference correction method was done in several 
steps. At first, type and initial location of the singularities representing each model were 
specified. Rules of thumb given in Appendix 9 were used to select type, location, and 
weighting factors for these singularities. A total of 11 singularities were selected for each 
model. A source and a sink were selected to represent fuselage blockage, two sources 
were selected to model the separation wake blockage effects. Seven line doublets, located 
along the 1/4 chord fine of the wing, were chosen to represent lifting effects. Weighting 
factors for the line doublets were selected to model an elliptic lift distribution for the wing. 
Figures 16a, 16b give the singularity representation of each semispan model for 0° angle of 
attack. The real-time coordinates of these singularities as a function of the pitch angle 
were computed using the known kinematics of a semispan model mounted on the image 
plane (for more detail see Eqs. (14. 19a), (14. 19b) in Appendix 14) . 

In the next step, measured lift force in combination with the Kutta-Joukowski formula 
was used to determine the strength of line doublets for the wing for each data point. 

The strength of the remaining singularities was computed using a least squares fit 
of the wall pressure measurements on 180 wall pressure ports that were arranged in six 
rows above the image plane. The least squares fit used wall pressure port rows 1,2, 3, 6, 7, 8 
depicted in Fig. 17a, 17b . For more detail on the least squares fit procedure see Section 2.4 . 

The standard deviation of the least squares fit of the wall signature was computed 
for each data point of Runs 154 and 219 (see Figs. 18a, 18b) . The standard deviation of 
the 8 % scale model was on the order of 0.002 in units of the dimensionless perturbation 
velocity. This agrees with the standard deviation of a wall signature obtained by Rueger et 
al ., [19] who reported a value of 0.005 in units of pressure coefficient, i.e. 0.0025 in units of 
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the dimensionless perturbation velocity. The standard deviation of the 14 % scale model 
was on the order of 0.003 to 0.006 in units of the dimensionless perturbation velocity. 

Figures 19a to 19f show the result of the least squares fit of the wall signature for the 
8 % semispan model at 19.82° uncorrected angle of attack. Figures 20a to 20f show the 
result of the least squares fit of the wall signature for the 14 % semispan model at 9.98° 
uncorrected angle of attack. Figures 21a, 21b depict the wall signature for both models at 
approximately 0.0° angle of attack at wall pressure port Row 6 . The large difference in 
solid volume blockage of both models can clearly be detected in the wall signature. The 
measured wall signature difference “u” depicted in Figs. (19a) to (21b), i.e. the velocity 
difference [U e + u t ] — U e in Fig. 9b, shows excellent agreement with its least squares fit. 

The present method (WICS), the two- variable method, and the classical method were 
used to compute wall interference corrections. Two-variable method results were provided 
by Mat Rueger of Boeing St. Louis. Classical corrections were provided by Alan Boone 
of NASA ARC who used NACA Rep. No. 995 (solid volume blockage), [20], R.A.E. Rep. 
No. 3400, [21] (separation wake blockage), and NACA TN 2454, [18] , to determine wall 
interference corrections. Mean wall interference corrections for each model were computed 
using flow field reference points located along the 3/4 chord fine of the wing. Corresponding 
results are compared below. 

As expected, wall interference corrections computed by WICS and the two-variable 
method show excellent agreement because both methods are based on potential flow theory 
and boundary flow measurements. Angle of attack corrections agree well in all three cases 
(see Figs. 22a, 22b) . The solid volume blockage factor contribution depicted in Fig. 23a 
agrees well for the 8 % scale model in all three cases. A comparison of the solid volume 
blockage factor contribution of the 14 % scale model depicted in Fig. 23b shows larger 
differences between classical corrections and WICS. This can be explained by the fact 
that the calculation of the solid volume blockage using the classical method (NACA Rep. 
No. 995) assumes that a wind tunnel of constant cross-section extends to far upstream 
and downstream of the semispan model. This assumption, however, cannot be justified 
anymore in the case of the 14 % scale model as the fuselage length is 16.69 [ft] and the 
length of the image plane is « 20.0 [ft] . The classical method will therefore overpredict 
the solid volume blockage effect for the 14 % scale model. The separation wake blockage 
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factor contribution for the 8 % scale model determined based on the classical method [21] 
is larger than the blockage factor computed using WICS or the two— variable method (see 
Figs. 23a) . This agrees with observations reported in the literature, [5], [19] . 

The 8 % and 14 % scale model have identical geometry. Therefore it is possible to 
compare the minimum of the blockage factor of both models by using a scale factor law 
(for more detail see Appendix 20). Results discussed in Appendix 20 demonstrate that 
blockage corrections computed with the present method (WICS) satisfy this scale factor 
law. Thus, wall pressure measurement accuracy and the solid volume description used by 
the present method are sufficiently accurate for computing blockage effects. 

In general, the ratio between measurements and unknowns of a least squares fit has 
to be large to take full advantage of its smoothing characteristics. In our application 
the number of unknowns of the least squares fit is two (see Eqs. (13a), (16a)). Figures 
23c, 23d compare the computed dynamic pressure correction for the 8 % scale model 
with 180 or 30 wall pressure ports used for the least squares fit of the wall signature. 
The differences in the computed corrections depicted in Figs. 23c, 23d are small. This 
demonstrates a key operational advantage of the Wall Signature Method : the calculation 
of the corrections is relatively insensitive to the number and location of the wall pressure 
ports (see also Table 1 in Chapter 1) . Comparison of the data scatter in the computed 
dynamic pressure correction depicted in Figs. 23c, 23d shows that an increase in the number 
of wall pressure measurements used in the least squares fit reduces the data scatter of the 
computed blockage corrections. 

The local dynamic pressure correction for the 8 % model at 19.82° angle of attack 
and for the 14 % model at 9.98° angle of attack are computed in the plane Y=0.0 [ft] (cut 
through test section parallel to side wall). The dynamic pressure correction for the 8 % 
model at 19.82° increases gradually as an observer moves from upstream to downstream 
of the model (see Fig. 24a) . The separation wake blockage effects dominate the dynamic 
pressure correction downstream of the model at an angle of attack of 19.82° . The contour 
lines are nearly parallel to the z-axis. We conclude that the dynamic pressure correction 
for the 8 % model is almost exclusively a function of the streamwise coordinate. 

The dynamic pressure correction for the 14 % model at 9.98° has a saddle point at the 
streamwise coordinate 120.0 [ft] (see Fig. 24b) . The correction decreases as an observer 
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moves downstream of the model. Therefore, it can be concluded that the solid volume 
blockage effects dominate the correction at an angle of attack of 9.98° . 

In both contour plots, it can be seen that the minimum of the dynamic pressure 
correction for a constant streamwise coordinate is about 4.0 [ft] above the image plane. 
This is caused by the fact that the tunnel width increases as an observer moves from the 
image plane surface to the tunnel centerline (see also Fig. 15). 

The local angle of attack corrections for 19.82° and 9.98° for each model are computed 
in the plane Y=0.0 [ft] (cut through test section parallel to side wall). The contour plots 
of the results axe depicted in Figs. 24c, 24d . The angle of attack correction variation along 
the 3/4 chord line of the 8 % model at 19.82°, i.e. from wing root to tip, is on the order 
of 0.05° . The aerodynamic twist is insignificant in this case. However, the angle of attack 
correction variation along the 3/4 chord line of the 14 % model at 9.98° is on the order of 
1.0° and cannot be ignored anymore. 

3.4 Application of the Method to the Ames Bipod 

In 1995 the Ames Bipod, a floor-mounted support system, was tested in the NASA 
Ames 12ft Pressure Wind Tunnel (PWT). Blockage corrections due to support system wall 
interference effects were computed using the Wall Signature Method (WICS) as outlined 
in Section 2.3.2 . 

The application of the Wall Signature Method to t ie Ames Bipod was done in several 
steps. At first, type, initial location, and weight of singularities representing the Ames 
Bipod were specified. Figure 25a shows the Ames Bipod geometry and lists type, initial 
location, and weight of these singularities. A total number of 27 singularities was selected 
for the support system. Singularities No. 1 to No. 18, i.e. 9 source/sink pairs, were chosen 
to represent solid volume blockage effects of the support system. Singularities No. 19 to 
No. 27 were chosen to model separation wake blockage effects. 

In a second step, strength values of these singularities were computed using a least 
squares fit of the wall pressure measurements. Wall pressures on Row 1 to Row 8 (see 
Figs. 17a, 17b) were measured during the support system calibration as a function of the 
total pressure in the settling chamber and the Mach number at a test section reference 
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point. Wall pressure measurements obtained during an empty tunnel calibration were also 
available. The least squares fit was applied to the difference between the wall signature 
measured during the support system calibration and the wall signature measured dining 
the empty tunnel calibration. Figures 25b, 25c show the result of the least squares fit on 
Row 2 and Row 4 . Measured wall signature difference and least squares fit show reasonable 
agreement. Larger differences between the measured wall signature and the least squares 
fit on Row 4 are caused by the fact that singularities cannot be placed too close to the floor 
for numerical reasons (see also Chapter 2.6). Finally, the blockage factor was computed on 
planes Y=0.0 [ft] (parallel to test section side wall), X=120.71 [ft] (parallel to test section 
inlet), and Z=0.0 [ft] (parallel to test section floor) . 

Figure 25d shows the blockage factor on plane Y=0.0 [ft] that was computed using 
the initial singularity location given in Fig. 25a . Figure 25e shows the blockage factor 
on plane Y=0.0 [ft] that was computed after the standard deviation of the least squares 
fit was minimized as a function of the singularity location (see Appendix 8). Comparing 
both contour plots it can be recognized that differences in the computed blockage factor 
are small, i.e. computed corrections are not very sensitive to the location of singularities 
as long as they are placed at the location of the support system. Comparing singularity 
locations before and after the minimization it can be noticed that sources No. 1 to No. 9 
and sources No. 19 to No. 27 have moved closer together. The minimization procedure 
has correctly deduced from the wall signature difference that the front post of the Ames 
Bipod has a significantly larger diameter than the pitch strut, i.e. most of the separtion 
wake of the Ames Bipod is caused by the front post. 

Figure 25f shows the blockage factor on plane X=120.71 [ft] . In this contour plot it 
can be seen that the blockage factor increases in any direction if an observer moves closer 
to the test section wall. A minimum of the blockage factor is located « 1.0 [ft] above the 
tunnel centerline. The blockage factor increases significantly if an observer moves closer 
to the test section floor. 

Figure 25g shows the blockage factor on plane Z=0.0 [ft] . As expected, the blockage 
factor increases if an observer moves from an upstream position to a downstream position. 
This observation is caused by the fact that the separtion wake blockage dominates blockage 
corrections downstream of the Ames Bipod. 
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CHAPTER 4 


CONCLUSION AND REMARKS 

A revised version of the Wall Signature Method was developed which allows the user to 
predict Mach number, dynamic pressure, and angle of attack correction due to wall inter- 
ference effects in three-dimensional subsonic wind tunnel testing of aircraft models. This 
improved formulation of the Wall Signature Method uses lift force, propulsion simulator 
thrust force, pitching moment, wall pressure measurements, empty tunnel calibration, sup- 
port system calibration, a simplified representation of the test article and support system 
in terms of singularities, and precalculated normalized solutions of the subsonic poten- 
tial equation expressed as normalized perturbation velocities to predict wall interference 
corrections at a model reference point in real-time. 

The method is applicable to complex wind tunnel and support system configurations. 
Wall interference corrections are found by taking the difference between a simplified rep- 
resentation of the wind tunnel flow field in terms of singularities and the corresponding 
free-air solution. Computational procedures were developed to predict solid body blockage, 
separation wake blockage, propulsion simulator blockage, and lift interference correction if 
a fullspan model or a semispan model is tested. 

The definition of normalized perturbation velocities used by the method has been 
improved to allow the user to take full advantage of the geometry modeling capabilities of 
a three-dimensional panel method code. A “global” least squares fit procedure of the wall 
signature was also introduced to improve the application of the Wall Signature Method 
in real-time. Optimal locations of singularites are defined by minimizing the standard 
deviation of the least squares fit of the wall signature. 

Experimental data obtained during tests of two semispam models mounted on an image 
plane in the NASA 12ft Pressure Wind Tunnel were applied to the modified Wall Signa- 
ture Method. Blockage and angle of attack corrections were computed for different angle 
of attack settings. In all cases, computed angle of attack corrections show good agree- 
ment with corresponding classical corrections. Computed blockage corrections are smaller 
than corresponding classical corrections. This result agrees with observations reported in 
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literature, [5], [19] . 

Experimental data recorded during the calibration of the Ames Bipod was also suc- 
cessfully applied to the method. 

Further experimental studies have to be conducted in the future to gain confidence in 
the method. 
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APPENDIX 1 


HIGHER ORDER CORRECTION FORMULAE 


The testing of an aircraft model at a high angle of attack in a three-dimensional 
subsonic wind tunnel can create large separation wake blockage effects. In this case, first 
order approximations of the Mach number and dynamic pressure correction as a function 
of the blockage factor e axe no longer sufficient. Second order approximations have to be 
derived. 

A second order approximation of the Mach number correction can be found using a 
Taylor series expansion. The Mach number is expressed as a function of a small change in 
the fluid velocity. Assuming that this velocity change AU is related to the flow velocity U 
and blockage factor e as : 

AU = t-U (1.1) 


we get the following Taylor series expansion of the Mach number : 


rrN „„ /rr . dM(U) UU) d 2 M(U ) (f f/) 2 

M(U + eU) = M(U) + , X ; + 


( 1 . 2 ) 


dU 1! ' d U 2 2! 

It is necessary to express the Mach number as a function of the fluid velocity U . Applying 
the energy equation we know : 

c p T t = c„ T + ij- (1.3) 

Combining the energy equation with the definition of the Mach number, i.e. 

M 2 = — — = -£=- (1.4) 

a 2 7 RT v ’ 

and with the relationship between specific heat at constant pressure, isentropic exponent, 
and Gas constant, i.e. c p (7 — 1) = 7 R , we get : 
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After some algebra we get for the first and second derivative : 
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Combining Eqs. (1.2), (1.6a), (1.6b) and rearranging terms we get : 
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Using the nomenclature introduced in Chapter 2.1 and assuming that M(U + e U) = M 0 
and M(U) = M e , we get the second order approximation : 
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Similar to the Mach number correction it is necessary to derive a second order ap- 
proximation of the dynamic pressure correction using a Taylor series expansion. 

Assuming that the velocity change AU is related to the flow velocity U and blockage 
factor e according to Eq. (1.1) we get the following Taylor series expansion for the dynamic 
pressure : 
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The dynamic pressure has to be expressed as a function of the fluid velocity U . Applying 
the energy equation (Eq. (1.3)), the Mach number definition (Eq. (1.4)), the ideal gas law, 

i.e. : 

P 


= RT 


( 1 . 10 ) 


and assuming isentropic flow, i.e. : 
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we get 


q(U) = 


[ R Pt 

l/(7-l) 

c p Tt — 

U 2 1 

V(7-l) JJ2 

t 

•o* 

3 

k 


o 

4* 

2 


( 1 . 12 ) 
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After some algebra the first and second derivative can be obtained as : 

d q(U ) q(U) 


d 2 q(U) 


d U 

q(u) 


u 


2 - M 2 (U) 


2-5 M 2 (U ) + (2 - 7 ) • M 4 (U) 


d U 2 U 2 

Combining Eqs. (1.9), (1.13a), (1.13b) and rearranging terms, we get : 
q(U + eU) - q(U) 


q(u) 


+ 

+ 


2 - M 2 (U) 
5 


M 2 {U) + y --M\U) 


(1.13a) 

(1.136) 


(1.14) 


Again, using the nomenclature introduced in Chapter 2.1 and assuming that q(U + e U) = 
goo, q(U ) = g e , and M(U ) = M e (t/), we get the approximation : 


?oo 


2 - M. 2 (£l) 


► 6 + 


4- w 2 (tr) + 


< 2 (1.15) 
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APPENDIX 2 


PANEL METHOD CODE MODIFICATIONS 
Introduction 

In general, it is necessary to use a three-dimensional panel method code to calculate 
normalized perturbation velocities caused by singularities if the Wall Signature Method is 
applied to a wind tunnel with a non-rect angular but constant cross-section. Unfortunately, 
commercially available panel method codes do not allow the calculation of the flow field 
of a singularity placed inside of a wind tunnel. However, a few modifications to a panel 
method code can be introduced which make it possible to solve this type of internal flow 
field problem and to compute perturbation velocities. 

A panel method code may be used to find the velocity potential of a given internal 
flow problem if the corresponding internal flow geometry, e.g. the geometry of the wind 
tunnel test section of constant cross-sectional area, is paneled. The specification of normal 
velocities at the wind tunnel inlet is also required (see Ashby et al. [15] for more detail on 
the application of a panel method code to internal flow problems; see Katz and Plotkin [22] 
for a more detailed description of three-dimensional panel method codes) . 

The boundary value problem of a singularity placed inside a wind tunnel flow field is 
depicted in Fig. 26 . The velocity potential has to fulfill the Laplace equation : 

v 2 [ 4 >oo + <f>t ] =0 (2-1) 


where <poo is the free-stream potential and <f> t is the wind tunnel potential due to the wind 
tunnel walls and the singularity. Zero normal flow has to be satisfied across the wind 
tunnel wall surface and so we get: 


d_ ' 

<9n 


&00 


+ 4>t 


= 0 


(2.2a) 


The normal velocity vector u-^ 


has to be specified by the user at the test section inlet: 


d _ ' 

dn 


fioo 


+ 4>t 


Uoo 


(2.2 b) 
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Based on the principle of superposition (see Fig. 26) it is possible to express the wind 
tunnel potential <j>t as the sum of the singularity potential <f> a and the wind tunnel wall 
potential <f> w : 

<}>t — <j> a + <f>w (2.3) 

Combining Eqs. (2.1), (2. 3) and knowing that the singularity potential itself fullfills Laplace’s 
equation we get: 

V 2 [<£oo + K ] = 0 (2.4) 

Combining Eqs. (2.2a), (2.3) and rearranging terms we get the boundary condition across 
the wind tunnel wall surface: 


8 _ 

dn 


too + tv 


d_ 

dn 


to 


(2.5a) 


Combining Eqs. (2. 2b), (2. 3) and rearranging terms we get the boundary 7 condition at the 
test section inlet: 

d [ ] d [ 1 

too 4" tw = Uqq t° (2.56) 

The velocity vector and the singularity potential <js c are known and so the boundary 
value problem given by Eqs. (2.4), (2.5a), (2. 5b) can be solved. This requires the modi- 
fication of the original panel method code such that the boundary conditions given by 
Eqs. (2. 2a), (2. 2b) are replaced by the boundary conditions given by Eqs. (2. 5a), (2. 5b) . 
Only the normal velocity component due to the singularity, i.e. 


has to be added to prescribed normal velocities at panel centroids. 

Finally, the flow field solution of a singularity placed inside a wind t unn el test section 
can be found by superimposing the solution of the boundary value problem given by 
Eqs. (2.4), (2. 5a), (2. 5b), i.e. + <f > w , computed using the modified panel method code 
with the known analytic solution of the singularity potential to ■ 

The modification of the boundary conditions based on Eqs. (2. 5a), (2. 5b) does not 
cause numerical difficulties if a singularity like a point source, point sink, or point doublet 
is selected. In these cases, the singularity potential <f>, 7 vanishes if the distance between 
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singularity and panel centroid is large and therefore normalized perturbation velocities can 
be compute with no restrictions. However, the boundary conditions have to be modified 
carefully if a horseshoe vortex or a semi-infinite line doublet is selected as a singularity. 
Panel method code PMARC, [15] models the internal flow geometry as a closed box and 
a horseshoe vortex or a semi-infinite line doublet will intersect panels defining the outflow 
conditions. Therefore, the flow field of a line doublet has to be computed by solving the 
flow field of a point doublet and this solution has to be made perfectly symmetric or anti- 
symmetric relative to the streamwise coordinate of the point doublet. Finally, a numerical 
integration is applied in the streamwise direction to obtain the flow field of a line doublet 
(for more detail see Appendix 4). 

Modifications of a panel method code for semispan test are similar to the procedures 
described above. The user only has to make sure that normal velocities induced on the 
image plane surface by a singularity located on the image plane surface, i.e. point sources 
and sinks representing the semispan model fuselage volume, are zero. 

Numerical Verification 

A slender Rankine body ( l/2r 0 « 10) is selected to verify the proposed panel method 
code modification if a source or sink is selected as a singularity (see Fig. 27a for a detailed 
desciption of the Rankine body geometry) . A point source of strength +1.0 [ft 2 ] is 
located at X=116.0 [ft], Y=2.0 [ft], Z=-3.0 [ft] and a point sink of strength —1.0 [ft 2 ] is 
located at X=126.0 [ft], Y=2.0 [ft], Z=-3.0 [ft] . The shape and surface pressure coefficient 
distribution of the corresponding Rankine body can be described in analytic form using 
polar coordinates, [23] . 

Assuming x\ — x 0 < x < x\ + (i2 — xi)/2, we obtain : 

siny?i/2 

r\ - r 0 : 

sin^>i 

X\ — x 

COS 02 1 = 

r l 

Cp = 1 - 4 sin 2 ~ + 3sin 4 ^ 

Assuming x\ + (x 2 — x\)/2 < x < x^ + x 0 , we obtain : 

sin <£> 2/2 


(2.6a) 


( 2 . 66 ) 

(2.6c) 



X — X2 
COS <p 2 = 

r 2 

(2.76) 

c p = 1—4 sin 2 ^ + 3 sin 4 ^ 

(2.7c) 

where the stagnation point distance x 0 and the halfbody radius r„ 
strength a of +1.0 are given as: 

for a point source 

H 

o 

II 

(2.8a) 

o 

II 

(2.86) 


In a first step it is necessary to panel the selected Rankine body using the known 
description of its shape in polar coordinates (see Fig. 27b) . The free-air solution of the 
surface pressure coefficient distribution of this geometry is calculated using panel method 
code PMARC, [15] . The result of this panel method code calculation is compared with 
the corresponding analytic solution given by Eqs. (2.6c), (2.7c) . Figure 28 shows excellent 
agreement between the panel method and corresponding analytic solution of the surface 
pressure coefficient distribution of a Rankine body. 

In the second step the Rankine body is placed inside a wind tunnel test section. The 
chosen wind tunnel geometry is similar to a test section configuration of the NASA Ames 
12ft Pressure Wind Tunnel which will be used for semispan model tests. Figure 29 shows 
the corresponding wind tunnel wall, support system and Rankine body p an eling Velocities 
are computed along eight rows on the wind tunnel wall using the original version of panel 
method code PMARC, [15] . 

Then, panel method code PMARC, [15] was modified to solve the boundary value 
problem given by Eqs. (2.4), (2.5a), (2.5b) . The velocity field due to a point source and 
point sink is used to represent the Rankine body (see also Appendix 3) . The Rankine body 
is replaced by corresponding point source and sink of strength +1.0 [ft 2 ] at X=116.0 [ft], 
Y=2.0 [ft], Z=-3.0 [ft] and -1.0 [ft 2 ] located at X=126.0 [ft], Y=2.0 [ft], Z=-3.0 [ft] and 
velocities are again computed on eight rows using the modified version of the panel method 
code. 

Figure 30 compares dimensionless velocities of both calculations on Row 4 . Velocities 
show excellent agreement to verify the proposed panel method code modifications if point 
sources or point sinks are used to represent volume effects of a test article. 
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A rectangular wing (s/ c = 3.56) with a NACA 0012 airfoil section at an angle of attack 
of 5° is selected to verify the proposed panel method code modification if a semi-infinite 
line doublet is selected as a singularity (see Fig. 31 for a detailed desciption of the wing 
geometry) . 

The rectangular wing is placed inside a wind tunnel test section. Figure 32 shows 
the corresponding wind tunnel wall, support system and wing paneling. Lift coefficient 
of the wing and velocities along eight rows on the wind tunnel wall are computed using 
the original version of panel method code PMARC, [15] . Panel method code PMARC 
was also used to compute the wall signature at angle of attack 0° . This wall signature 
is due to the thickness of the wing and was therefore subtracted from the wall signature 
calculated for 5° to obtain the wall signature due to lift only. 

Then, panel method code PMARC, [15] was modified to solve the boundary value 
problem given by Eqs. (2.4), (2. 5a), (2. 5b) . The velocity field of four line doublets is used 
to represent lifting effects of the wing at 5° angle of attack (see also Appendix 4) . The 
location of these line doublets is depicted in Fig. 33 . The computed lift coefficient of 
the wing placed inside the wind tunnel test section (see Fig. 32) has to be related to the 
strength of the semi-infinite line doublets. Using the definition of the lift coefficient and 
the Kutta/Joukowski formula we get for the lift force : 

U 2 

L = cl • p • — — — • s-c — p-U-T-s (2.9) 


So we get for the circulation T : 

r 



( 2 . 10 ) 


The line doublet strength a of each of the four line doublets is then given as (see Ap- 
pendix 4) : 


er 




(2.11a) 


or 



(2.116) 


Knowing that cl = 0.3556, c = 1.5 [/<], and s/4 = 1.3334 [/<] for the selected wing 
we get for the singularity strength per unit velocity of each of the four line doublets 
a/U = 0.3556 [ft 2 ] . 
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The wing can now be replaced by four line doublets of strength 0.3556 [ft 2 ] and 
velocities are again computed on eight rows using the modified version of the panel method 
code. 

Figure 34 compares dimensionless velocities on the wind tunnel wall of both calcula- 
tions at Row 8. Velocities due to lift effects show excellent agreement to verify the proposed 
panel method code modifications if semi-infinite line doublets Eire used to represent lift of 
a wing. 


Summary 

A three-dimensional panel method code was modified to CEilculate the flow field of 
a singularity placed inside a wind tunnel test section. The modification was verified by 
replacing a Rankine body by a point source and point sink and by replacing a rectangular 
wing at 5° angle of attack by four line doublets. Corresponding flow field solutions compare 
favorably in both cases. The panel method code modification will work with no restrictions 
if a point source or point sink is selected as a singularity. However, boundary conditions 
have to be modified carefully if a horseshoe vortex or semi-infinite fine doublet is selected as 
this type of singularity intersects outflow panels of the internal flow geometry and therefore 
requires a numerical integration of the corresponding point doublet solution. 
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APPENDIX 3 


POINT SOURCE VELOCITY VECTOR 


The velocity vector of a three-dimensional point source is required if a panel method 
code is modified to compute the flow field of a source or sink placed inside a wind tunnel. 
The potential of a point source at location (15,2/5,25) is given as, [22] : 


<f>s(x,y,z) = 


— cr 


4 • 7r • [ [x - X 5] 2 + [y - ys ] 2 + [z — Z5] 2 ] 1/2 
The velocity components of this point source axe : 

d(j>s cr x — xs 


u(x,y,z) = 

v(x,y,z) = 
w(x,y,z) = 


dx 

d<h 

dy 

d<f> s 


4 v [ [* - X5] 2 + [y - ys ] 2 + [z - zs] 2 ] 
cr y-ys 


3/2 


4 n [ [x - X5] 2 + [y - ys] 2 + [z - zs] 2 ] 
o z - zs 


(3.1) 


(3.2a) 


575 ( 3 - 2t o 

(3.2c) 


8 * 4 1 [ {* - *s] 2 + [» - ys] 2 + [z - 2 S ] 2 ] 3 ' 2 

The singularity strength cr has the unit [ft 3 /sec] or [m 3 /sec] . The perturbation 
velocity field given in Eqs. (3. 2a), (3. 2b), (3. 2c) could be used to change a panel method 
code as outlined in Appendix 2 . 
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APPENDIX 4 


LINE DOUBLET VELOCITY VECTOR 


The velocity vector of a three-dimensional semi-infinite line doublet is required if a 
panel method code is modified to compute the flow field of a line doublet placed inside a 
wind tunnel. 

In general, a semi-infinite line doublet can be considered as an elementary horseshoe 
vortex ([22] , see also Fig. 35) . The potential of a line doublet cam be obtained by integrating 
the solution for a point doublet in the x-direction, [24] . A point doublet does not have a 
radial symmetry as in case of a point source and therefore the line doublet potential will 
be a function of the orientation of the point doublet. For lift force acting in the positive z- 
direction (see Fig. 35) it is necessary that point doublets point in the negative z-direction. 
The potential of a finite length line doublet with starting point at (*i , yi , z\ ) and end point 
at ( XN,yi,zi ) is then given as, [22] : 




a 

4 7T 



[z - zi] 

[x-£] 2 + [y-yi] 2 + [z-zx] 2 ] 3/2 


(4.1) 


or 


<f>(x,y,z) ' = 

A = 


—a 


z — Z\ 


4 it [y- yi] 2 + [z- zi] 2 
x — Xjy 


[[x-xjv] 2 + [y-yi] 2 + [z-zi] 2 ] 

X — X\ 


1/2 


(4.2a) 


(4.26) 


[[x-xi] 2 + [y-yi] 2 + [z-zi] 2 ] 1/2 

Taking the limit xn —*■ oo we get the potential for the semi-i nfini te fine doublet : 


<j) L {x,y,z ) = lim <f>(x,y,z ) 

X ft — ►OO 


cr 

4 7T 


Z — Z i 

[y-yi ] 2 + [z-zi ] 2 


• B 


(4.3a) 


B = 1 + X —^ m (4.36) 

[[z-n] 2 + [y — yi} 2 + [z-*i] 2 ] ' 

This equation agrees with equation (5-35) in Ref. [24], if the starting point of the line 
doublet is located at (xi = 0,yi = 0,zi = 0) . The velocity components can now be 
obtained by talcing derivatives of the velocity potential <f>L • 
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In the x-direction we get : 


u(x,y,z) 


d<j>L 

dx 


In the y-direction we get : 


cr z — z\ 

4 * [[ x - xi ] 2 + [ y - yi ] 2 + [z-zi ] 2 ] 3/2 


(4.4) 


v(x,y,z) 


9^l 

dy 


—a 
4 7T 


[y - yi] • \z_ - z\] r 

[y-yi ] 2 + [z-z \] 2 L 


(4.5 a) 


c = 


B 


+ 


X — X\ 


[» - »') 2 + I z - 2 >) 2 [ [* - 111 2 + bf - w] 2 + [z - *,] 2 ] 3/2 

In the z-direction we get : 


(4.56) 


w(x,y,z) 


d < f> L 

dz 


cr 

4 TC 


1 

[y-yi ] 2 + [z - *i ] 2 


■ D 


(4.6a) 


d = ty-yi 3 2 - [*-*i] 2 . R 

[y-yi ] 2 + [z-zi ] 2 

_ [x - Xi] • [z - Zi ] 2 (4.66) 

[[x-xi] 2 + [y-yi ] 2 + [z-zi ] 2 ] 3/2 

The caJculation of the flow field of a semi-infinite line doublet located inside a wind 
tunnel of constant cross-sectional area is difficult if a modified version of panel method 
code PMARC is used (see Appendix 2) . The line doublet intersects the exit plane of the 
paneled wind tunnel geometry and can cause convergence and accuracy problems. 

Fortunately, it is possible to calculate the flow field of a fine doublet in a wind t unne l 
of constant cross-sectional area by integrating the flow field of the panel method code 
solution of a point doublet in the stream wise direction. The line doublet flow field solu- 
tion is obtained by shifting and superimposing corresponding point doublet solutions (see 
Fig. 36a) . 

In general, the potential of a semi-infinite line doublet (Eq. (4.1) ; xjv >> ®i) c an be 
approximated using numerical integration. We get then : 


N 

< t > L { x , y > z ) ss ^ 2 < t > pD ( x , y , z ; x* ) • Ax 
Jfc=i 


(4.7a) 


< f > PD ( x > y,z ; x fc ) 


cr z — z \ 

47r [[x-x fc ]2 + [y-yi ] 2 + [z- 2l ]2] 3/2 


(4.76) 
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Sfc = *l + (fc — 1)* Ax 


XN » Xl 


where Ax is the step size of the integration. The perturbation velocities of the line doublet 
can now be approximated as : 

N 


u(x,y,z) 

v(x,y,z) 

w(x,y,z) 


d 


Ax ' foT <I>pd(x, y, z ; x* ) 


jt=i 

N 


Ax ' S IT ^ PD ^ X ’ y ’ z; Xk ) 

jt=i oy 


N 


d 


Ax ' S ^ pD ( x » y» z ; x * ) 


it=i 


(4.8a) 

(4.86) 

(4.8c) 


where 


— <j>pD(x, y, z; xjt ) 


(-3) • [x - x fc ] • [z - zi] 


4 7T 


dy 


<f>PD(x,y, z; Xjt ) = 


— <^>pd(x, y, z; xjt ) 


[[x-xjt] 2 + [y — yi] 2 + [z-zi ] 2 ] 5/2 

£ . (—3) • [y — yi] • [z — z\\ 

4r [[ I -**] 2 + [y- W P + [z-z,] 2 ] 5 ' 2 
cr [x - xjt] 2 + [y - yi] 2 - 2 • [z - z^ 2 

,5/2 


(4.9a) 

(4.96) 

(4.9c) 


4 x [ [x - Xjt ] 2 + [y-yi ] 2 + [z-zi ] 2 ]' 

For a wind tunnel of constant cross-sectional area a streamwise shift of the flow field 
solution of a point doublet located at (xi,yi,zi) is possible (see Fig. 36b) . Therefore 
Eqs. (4.8a), (4.8b), (4.8c) only require the calculation of <j>PD{x,y,z ; ii ) as the following 
relationships apply : 

d d 

— (f>PD(x,y,z-, xjt ) = 4 >pd{x — [£ - l]Ax, y,z; xi ) ( 4 . 10 a) 

d d 

— <f> PD (x,y,z; Xjt ) = -7^ 4 >pd(x - [fc - l]Ax,y,z; xi ) ( 4 . 106 ) 

d d 

<f>PD(x,y,z; x k ) = -qJ <f>PD(x - [fc - l]Ax, y,z; x x ) ( 4 . 10 c) 

The numerical integration defined by Eqs. (4.8a), (4.8b), (4.8c) has to be done care- 
fully. Integration error due to the fact that the numerical solution of a point doublet flow 
field is not perfectly antisymmetric (see Eq. (4.9a)) or symmetric (see Eqs. (4. 9b), (4.9c)) 
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relative to its x-coordinate x* has to be avoided. This can be done by considering, e.g., 
the numerical solution upstream (x < ifc) of the point doublet as exact and imposing the 
following conditions on the downstream part (x > xjt) of the perturbation velocities : 


u(x,y,z) = 

- u(x k -[x-x k ],y,z) ; 

X > Xjfc 

(4.11a) 

v(x,y,z) = 

v(x k - [x - Xjfc],y. z) ; 

x > x k 

(4.116) 

w(x,y,z) = 

w(x k - [x - xjfe],y,z) ; 

x > x k 

(4.11c) 


Unfortunately, the point doublet and therefore also the semi-infinite fine doublet does 
not have a radial symmetry. However, the velocity vector of a point or line doublet not 
pointing in the negative z-axis direction can easily be found by applying a coordinate 
system rotation. A rotation angle r related to the directional property of the fine doublet 
is introduced as depicted in Fig. 37 . 

The relationship of coordinates and velocities between the line doublet fixed coordinate 
system (x,y,z) and reference coordinate system (x\y , ,z') is depicted in Fig. 37 . The 
rotation angle direction is defined such that the lift force caused by the semi-infinite fine 
doublet points in the positive t/-axis if the line doublet is rotated by +90° . In this 
case we get for the transformation of coordinates from the reference coordinate sytem, i.e. 
(x / ,y / ,z / ), to coordinates in the line doublet fixed coordinate system (x,y,z) : 

x = x' (4.12a) 

y = y' • cosr - z'-sinr (4.126) 

z = y'-sinr + z'-cosr (4.12c) 

The velocities components (u,u,ty) axe computed i:i the line doublet fixed coordinate 
system. Finally it is necessary to back-transform these velocity components to the refer- 
ence coordinate system. We get then : 


u — u 


v' = v ■ cos r + w • sin r 


(4.13a) 

(4.136) 
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— t; • sin r + w • cos r 


w 


(4.13c) 


The singularity strength cr has the unit [ft 3 /sec] or [m 3 /sec] . The perturbation 
velocity field given in Eqs. (4. 8a), (4.8b), (4.8c) could be used to change a panel method 
code as outlined in Appendix 2 . 
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APPENDIX 5 


CALCULATION OF THE LINE DOUBLET STRENGTH 

In general, it is possible to compute the strength of line doublets representing lifting 
effects of an aircraft by using lift and pitching moment measurements recorded during a 
wind tunnel test (see Ulbrich and Sfezn/e,[7] ). 

Assuming that the lifting surfaces of the wing and tail of an aircraft model axe dis- 
cretized by using equally spaced line doublets along the 1/4-chord line of the wing and 
tail we get for the total lift and pitching moment (see also Fig. 38) : 


Tl u 

L = £l„( 

>) + f>y) 

(5.1) 

t = l 

j=l 


fl w r 

n t 


P — ^ ^ L w {%) • [ x mr x t u (i) ] 

+ Lt (fi '{ Xmr ~ **(*) i 

(5.2) 

t=i L 

i= i L J 



where L w (i ) and L t (j ) are the lift contributions of the line doublets of the wing and tail. 

Similar to Eqs. (2.9) and (2.11), the Kutta/Joukowski formula may be used to connect 
a lift force to a line doublet strength. The appli cation of the Kutta/Joukowski formula in 
a wind tunnel requires the calculation of the product p^ • U a 0 where poo is the free-stream 
density and I/oo is the free-stream velocity at the location of the aircraft model. A first 
order approximation of the product poo • Uoa may be obtained by applying wall interference 
corrections to the measured reference density p re f and reference velocity U re f at a model 
reference station (see Eqs. (5.78) and (5.82) in Ref. [10]) . We then get : 


Poo * Uoo ~ 

Pref ‘ U re f * 1 + 6 • ( 1 M ^ e f ) 


(5.3a) 

The Kutta/Joukowski formula may now be written as : 



— Poo Uoq (0 

~ Pref * Uref ’ [ ^ C ° ( 1 — ) 

• <(*) 

(5.36) 

Lt(j) ~ Poo * Uqo ° (i) 

~ Pref * Uref * £ 1 + £ * ( 1 — ) 

■ <(j) 

(5.3c) 


Introducing singularity weights 1 ^( 1 ), wt(j), and discrete line doublet span A s w , A s t 
on the wing and tail we get for the singularity strength : 


0 -* (j) = T w -As w -W w (i) (5.4a) 


57 



r* • a s t • wt(j) 


(5.4b) 


<( 3 ) = 


The circulation and T t axe the only unknowns in Eqs. (5.1) and (5.2) if line doublets 
weights are specified by the user. The weights may represent, e.g., elliptic lift distribution 
along the wing and tail of an aircraft model. Experimental or CFD solutions of the true 
lift distribution may also be used to determine weights (see also NACA Report No. 921, 
[25]). Combining Eqs. (5.1) to (5.4b) we get : 

L 


Pref • U ref • [ 1 + e • ( 1 -M r 2 e/ ) 
p 

Pref ' Uref * 1 + C • ( 1 — M ) 


r w • Q£i + r t • «2 


— Tv-fa -f r t • fa 


where 


Tl w 


qi = a s w • y; w w (j) 


«•= 1 


a 2 = A s t y ^wt(j) 
j = i 


Ti w 

fa = A s w • 


f=i L 
n t r 


w w (t) • [ x mr x w (i) ] 
wt(j) ■ [ %mr x tU ) 1 


(5.5a) 

(5.56) 

(5.6a) 

(5.66) 

(5.7a) 

(5.76) 


fa = As t • ^ 

3= 1 L 

Equations (5.5a) and (5.5b) are a 2 x 2 linear system of equations which can be solved 
easily for T*, and T t : 

L • fa — P • a 2 


= 


r« = 


Pref • U re f - [ 1 + e • ( 1 - M 2 e/ ) ] 0‘1-fa - 0‘2-fa 

(-L) -fa + P - ai 


(5.8a) 


(5.86) 


Pref'U ref • [l + e-( 1 -M 2 e/ )] <* t -fa - a 2 -fa 

Finally, the singularity strength is obtained by applying Eqs. (5.4a) and (5.4b) . 

The line doublet strength has to be compatible with non-dimensionalized perturba- 
tion velocities used in the panel method code computation of the wind tunnel and wall 
interference flow field (see Eq. (25) and Appendix 2). Correct units for velocity, density, 
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force, and pitching moment have to be selected if Eqs. (5.4a), (5.4b) are used. Table 4 lists 
units : 


Table 4: Units for Singularity Strength Calculation 


Quantity 

SI unit 

non-SI unit 

Velocity 

Density 

Force 

Moment 

m / sec 
kg j m 3 

N = kg ■ m / sec 2 
N • m 

ft j sec 
Slug / ft 3 

Ibf = Slug • ft / sec 2 
ft ■ Ibf 


The singularity strength as defined in Eqs. (5.3b), (5.3c) has the unit [m 3 /sec] or 
[ft 3 ! sec] . Similar to Eq. (25), it is necessary to divide the strength by the reference 
velocity U re f • Therefore we get for the normalized strength in units [m 2 ] or [ft 2 ] : 

MO = (5.9a) 

a,(j) = («•«) 

Uref 

Semispan model tests are often conducted using only the wing and fuselage of an 
aircraft model. In this case the circulation T t is zero and only the lift force measurement 
will be used to calculate the strength of line doublets representing lifting effects of the 
wing. Using Eqs. (5.5a) we get for : 


r* 


L 


Pref ' U re f ■ 


! + «•(! ~K. t ) 


* O'! 


(5.10) 


In general, the application of Eqs. (5.8a), (5.8b), (5.10) to the Wall Signature Method 
requires at least a one step iteration as blockage factor t can only be estimated after 
a successful least squares fit of the wall signature. Zero may be selected as an initial 
approximation of blockage factor e in Eqs. (5.8a), (5.8b), and (5.10). 
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APPENDIX 6 


TRI-LINE AR INTERPOLATION 


In general, tri-lineax interpolation can be used to find the value of a scalar function 
f(x p ,y p ,z p ) using known values f(x{ , yj , z k ),f(x i+ i , yj , z k ), ■■■ , at corresponding eight 
grid cell corner points (see Fig. 39) . 

Tri-linear interpolation requires several steps. At first we have to find grid cell indices 
(i,j,k) such that x t - < x p < Xi+i ; yj < y p < yj + 1 ; z k < z p < . For a constant grid 

cell size Ax, Ay, A z, indices (i,j, k ) can be computed as : 


+ 1 

(6.1a) 

int[ yp ~ yi ] +1 

(6.16) 

+ 1 

(6.1c) 


k = 

where xi, yi, z\ are the starting coordinates of the entire grid. 

In the next step, weighting factors of each grid cell comer point have to be found. 
The fine connecting a corner point with its opposite comer point is assumed to be a 
spacial diagonal of the grid cell. Then, the weighting factor of a grid cell comer point is 
computed by dividing the volume of the rectangular prism defined by point (x p ,y p ,z p ) and 
the opposite corner point by the total volume of the grid cell. For example, the weighting 
factor of grid cell corner point (x t +i , yj ,z k ) is : 

K',i+i,Jfc+i 


w i+i,j,k = 


V c 


cell 


(6.2a) 


where 


K',j+i,jt+i — ABS 


[x p - x t -] • [y p - y j+ 1 ] - [z p - z k+1 ] ] (6.26) 

V ce ii = ABS [x l+ i - x t ] • [y J+ 1 - yj] • [z k +i - z k ] ] = Ax • Ay • A z (6.2 c) 

Finally, the functional value f(x p , y P ,z p ) can be interpolated as : 


t+i j + l fc+i 

f( x p>yp> z p) = SEE a Xot j VP > ^7 ) 

a=t P 7 = fc 


(6.3) 
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Table 5 lists grid cell corner points and corresponding opposite comer points. 


Table 5: Grid Cell Indices 


No. 

Corner Point 

Opposite Corner Point 

1 

h j j k 

i + l t j + i,k + l 

2 


i+l,j,k + 1 

3 

i + l,j,k 

i,j + l,k + l 

4 

* + 1 j J + 1 1 

+ 1 

5 

i,j,k + 1 

i + 1, J + 1, k 

6 

i,j + l,fc + 1 

*' + 1 J,k 

7 

i + l,j,k + 1 

i>j "b k 

8 

* + l,j + l,fc + l 

i,j,k 


Tri-linear interpolation is applied to calculate normalized perturbation velocities of a 
singularity located at a point (x p ,y p ,z p ) assuming that perturbation velocities axe known 
for singularities located on eight comer points of a singularity grid cell. 
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APPENDIX 7 


PARABOLIC INTERPOLATION 


A parabola is used to interpolate normalized perturbation velocities y(M) as a function 
of the real-time Mach number M . The paxabola has the following form : 


y{M) = a ■ M 2 + b ■ M + c 


(7.1) 


Three coefficients a, b, c of the parabola have to be calculated. Therefore it is necessary 
to compute the perturbation velocity database for three discrete Mach numbers. 

Three discrete Mach numbers M\ — 0.0, M 2 = 0.3, M 3 = 0.6 axe selected and corre- 
sponding normalized perturbation velocities y(Mi), y(M 2 ), y{Mz) axe known. Then we 
get for the coefficients of the paxabola : 


a 


1 

M3 — M2 


y(M 3 ) - y(Mi) 

M3 — Mi 


y(M 2 ) - y(Mi) 
M 2 — Mi 


(7.2a) 


6 


y(M 2 ) - y(Mi) 
M 2 — Mi 


a • (Mi + M 2 ) 


(7.26) 


c = y(Mi) — a • M 2 — b-M\ 


(7.2c) 


Coefficients a, 6 , c have to be computed for ni • (n 2 + ns) parabolas in real-time if, 
e.g., a test article and support system axe represented by ni singularities, wall signatures 
axe measured at a total number of n 2 wall pressure ports, and a reference point grid of 
points is selected. 
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APPENDIX 8 


OPTIMIZATION OF THE SINGULARITY LOCATION 

The wall signature method calculates wall interference corrections based on a singu- 
larity representation of the test article and the support system. The wind tunnel test 
engineer has to specify the initial location of these singularities using simple rules. For 
example, sources and sinks representing the fuselage volume should be located along the 
fuselage axis such that they represent a Rankine body approximation of fuselage volume 
effects, line doublets should be distributed along the 1/4-chord line of lifting surfaces, and 
sources related to separation wake effects should be located at points where separation is 
likely to occur during a test. 

The initial singularity representation is not unique as the test engineer has to specify 
the singularity location. However, it is possible to make the singularity location unique by 
minimizing the standard deviation of the least squares fit as a function of the singularity 
location. 

An efficient optimization can only be achieved if the total number of independent 
variables of the minimization, i.e. the coordinates of each singularity, is reduced to a 
reasonable limit. It is also necessary to restrict the direction of change of the singularity 
coordinates. 

During the support system calibration three groups of singularities exist, i.e., sources of 
the support system volume, sinks of the support system volume, and sources of the support 
system wake. Therefore it is possible to reduce the number of independent variables 
down to three, if the relative distances of the singularities of each group is kept constant 
during the optimization. The direction of change of the optimization is defined by the 
line connecting the first source and the first sink of singularities representing the support 
system volume. 

Similarly, during the wind tunnel test of a test article five groups of singularities 
exist, i.e., sources of the fuselage volume, sinks of the fuselage volume, and sources of the 
wing separation wake, line doublets of the wing, and fine doublets of the tail. Therefore 
the number of independent variables can be reduced to five, if the relative distances of 
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the singularities of each group is kept constant during the optimization. In this case the 
direction of change of the optimization is defined by the fine connecting the first source 
and the first sink of singularities representing the fuselage volume. Assuming that the 
1/4-chord line of the wing and tail Eire ideal locations for line doublets of the wing an d 
tail it is even possible to reduce the number of independent variables to three. 

An optimization algorithm based on the Method of Steepest Decent has been included 
in the real-time software package of WICS (see also Ref. [26]) . 
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APPENDIX 9 


SELECTION OF THE SINGULARITY / REFERENCE POINT LOCATION 

The Wall Signature Method uses a singularity representation of the wind tunnel model 
to predict wall interference corrections. The singularity representation of the model is 
uniquely defined if TYPE, LOCATION, and STRENGTH of each singularity are known. 

The STRENGTH of each singularity is derived from real-time measurements of wall 
signature, lift force, and pitching moment. If a propulsion simulator is used during a test 
it is also necessary to measure the propulsion simulator thrust and the propeller disk area. 

The TYPE of each singularity is chosen by the wind tunnel test engineer. A singularity 
TYPE should be selected such that it represents a blockage or lifting effect of the wind 
tunnel model. Sources and sinks are used to represent volume and wake blockage effects 
of the wind tunnel model. Line doublets are used to represent lifting effects of the wind 
tunnel model. A sink is used to model blockage effects of a propulsion simulator in wind 
tunnel testing. 

The LOCATION of each singularity has to be specified by wind tunnel test engineer. 
LOCATIONS should be selected based on geometry of wind tunnel model. All singularity 
coordinates have to be provided in tunnel coordinates. 

The following empirical rules will help the test engineer to make reasonable selections 
of the singularity TYPE and LOCATION : 

(1) FULLSPAN WIND TUNNEL MODEL : 

(1.1) Fuselage Volume : A source of weighting factor “+1.0” has to be placed on 
the fuselage axis approximately one mean fuselage radius downstream of the nose of the 
fuselage. A source of weighting factor “-1.0” , i.e. a sink, has to be placed on the fuselage 
axis approximately one mean fuselage radius upstream of the tail end of the fuselage. 

(1.2) Wake Separation : Sources with positive weighting factors have to be placed on 
the wind tunnel model at locations where flow separation occurs. If a wind tunnel model 
is tested, e.g., in landing configuration sources of equal strength should be placed at the 
location of the wing flaps. 
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(1.3) Wing : The wing span is divided into equal size wing span increments. One line 
doublet is assigned to each wing span increment. Line doublet starting points are placed 
where the 1/4 chord line of the wing and the middle of each wing span increment intersect. 

(1.4) Tail : The tail span is divided into equal size tail span increments. One line 
doublet is assigned to each tail span increment. Line doublet starting points are placed 
where the 1/4 chord line of the tail and the middle of each tail span increment intersect. 

(1.5) Propulsion Simulator : A source of weighting factor “-1.0” , i.e. a sink, is placed 
at the center of the propeller if a turboprop engine is simulated; a sink is placed halfway 
between the compressor and turbine if a turbojet or turbofan engine is simulated. 

(2) SEMISPAN WIND TUNNEL MODEL : 

(2.1) Fuselage Volume : A source of weighting factor “+1.0” has to be placed on 
the fuselage axis approximately one mean fuselage radius downstream of the nose of the 
fuselage. A source of weighting factor “-1.0” , i.e. a sink, has to be placed on the fuselage 
axis approximately one mean fuselage radius upstream of the tail end of the fuselage. 
The z-coordinate of the source and sink of the semispan model fuselage does not have to 
be specified. The Wall Interference Correction System (WICS) of the NASA Ames 12ft 
Pressure Wind Tunnel (PWT) implicitly assumes that their z-coordinate is identical with 
the z-coordinate of the image plane surface. 

(2.2) Wake Separation : Sources with positive weighting factors have to be placed on 
the wind tunnel model at locations where flow separation occurs. If a wind tunnel model 
is tested, e.g., in landing configuration sources of equal strength should be placed at the 
location of the wing flaps. 

(2.3) Wing : The wing semispan is divided into equal size wing span increments. One 
line doublet is assigned to each wing span increment. Line doublet starting points are 
placed where the 1/4 chord line of the semispan wing and the middle of each wing span 
increment intersect. 

(2.4) Propulsion Simulator : A source of weighting factor “-1.0” , i.e. a sink, is placed 
at the center of the propeller if a turboprop engine is simulated; a sink is placed halfway 
between the compressor and turbine if a turbojet or turbofan engine is simulated. 
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(3) REFERENCE POINTS : 

WICS computes mean wall interference corrections for sets of reference points. The 
test engineer has to specify these sets of reference points in the tunnel coordinate system. 
WICS allows the test engineer to specify up to 10 independent sets of reference points. 
Sets of reference points can be specified along the fuselage axis and along the 3/4 chord 
line of the wing or tail. 

A simple Expert System should be designed in the future that assists the test engineer 
in the selection of the singularity location and weighting factors. 
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APPENDIX 10 


BLOCKAGE EFFECT OF A POWERED WIND TUNNEL MODEL 


The correct simulation of flow interference effects between the aircraft fuselage, the 
wing, and the engine cowling during a wind tunnel test requires the installation of a propul- 
sion simulator in a wind tunnel model. The blockage correction caused by a propulsion 
simulator has to be estimated if a significant amount of thrust is produced. 

In general, a propulsion simulator is operated as a small propeller. A point doublet 
pointing in the streamwise direction may be used to estimate propulsion simulator blockage 
effects during a wind tunnel test, [22] . A more accurate modeling of blockage effects may 
be obtained if a semi-infinite line doublet pointing in the streamwise direction is used, [27] . 
The velocity potential of a semi-infinite line doublet with a starting point at (xi , y \ , z \ ) is 
given as, [22] : 




— <T 

4 7 r 



[*-a <jj 

>-{P + fo-mP + [z-*iP] 3/! 


( 10 . 1 ) 


It can be shown that the velocity potential defined in Eq. (10.1) is identical with the 
velocity potential of a sink located at (®j, yi,zi) (see Eq. 3.1) : 


<f>LD(x>y,z) = - <j>s{x,y,z ) = 


< 7 

4 7 r 


1 

[l*-*iP + fe-mP + [z-*iP ]‘ /2 


( 10 . 2 ) 


Location and strength of this sink have to be specified. It is reasonable to place the sink 
at the center of the propeller or compressor disk of the propulsion simulator. The sink 
strength a may be related to the propulsion simulator thrust Tj>, propeller disk area 5, 
free-stream velocity , and free-stream density p^ by applying an approximate solution 
of the flow around the ideal propeller given by Koning , [6] . In wind tunnel testing, free- 
stream velocity and density have to be estimated by applying blockage corrections to the 
measured reference velocity t/ re / and reference density p re f, i-e. Uoo ~ U re f (1 + 0 and 
Poo ~ Prej (1 — e ^ref ) ( see Eqs. (5.78) and (5.82) in Ref. [10]) . Then, we get for the 
propulsion simulator sink strength : 


= 5 



(1 + e ) 2 + 


2 ■ Tp 


w(i-<«y-s 


- U„j •(! + £) 


(10.3) 
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Equations (10.3), i.e. Koning’s formula, may also be obtained by considering blockage 
effects of a lifting rotor in subsonic wind tunnel testing, [27] . 

Effects of a propulsion simulator on wind tunnel blockage effects can easily be in- 
cluded in the least squares fit of the wall signature required for the application of the 
Wall Signature Method. The strength of the propulsion simulator sink is known from the 
thrust measurement. It is only necessary to subtract the corresponding wall signature 
contribution from the total wall signature (see also Eqs. (13g) in Chapter 2) . 

The sink strength <r* as defined in Eq. (10.3) has the unit [m 3 /sec] or [ft 3 /sec] . It 
has to be divided by I/ re / if normalized perturbation velocities as defined in Chapter 2 are 
used for the least squares fit of the wall signature. Finally we get : 


= S- 


2 • Tp 


(1 + e) 2 + — 2 j- — - (1+e) 

Pref ■U? er S-(l-e M 2 e/ ) 


(10.4) 


The application of Eq. (10.4) to the Wall Signature Method requires at least a one 
step iteration as blockage factor e can only be estimated after a successful least squares fit 
of the wall signature (see also Appendix 5). 
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APPENDIX 11 


QUALITY CHECK OF THE LEAST SQUARES FIT 

The Wall Signature Method uses a linear least squares fit of the difference between 
wall signature and corresponding wall pressure port calibration to predict blockage effects 
of a test article. In practical applications the facility hardware cannot guarantee that good 
pressure measurements are recorded on all wall pressure ports at all times. However, not 
all pressure measurements are required to compute a blockage correction of the test article. 
Therefore an efficient implementation of a wall interference correction system based on the 
Wall Signature Method has to check measurements on each wall pressure port. 

Studies using experimental data have shown that a four— step quality check of the 
least squares fit is sufficient to identify unacceptable wall signature measurements without 
rejecting too many data points. Wall signatures obtained during the empty tunnel or 
support system calibration have to be inspected separately. They are used to remove 
orifice error, wall divergence, and wall boundary layer displacement effects. 

The proposed quality check of the difference between the real-time and calibrated 
wall signature is done as follows (see also Fig. 40) : 

CHECK 1 : WALL SIGNATURE DIFFERENCE MAGNITUDE 

The absolute value of the wall signature difference is computed for each wall pressure 
port. A port is rejected, i.e. its port flag is set to zero, if the wall signature difference is 
larger than a specified upper bound. This bound is a function of the wind tunnel facility 
configuration and data aquisition hardware. An upper bound of 0.1 has been selected for 
the NASA Ames 12ft Pressure Wind Tunnel. 

CHECK 2 : OUTLIER IDENTIFICATION 

The absolute value of the difference between the wall signature difference and its least 
squares fit is computed for each wall pressure port. A port is rejected, i.e. its port flag is 
set to zero, if this difference is larger than three times the standard deviation of the least 
squares fit. 

CHECK 3 : STANDARD DEVIATION OF EACH ROW 

The standard deviation of each wall pressure port row is computed. A wall pressure 
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port row is rejected, i.e. its wall pressure port flags axe set to zero, if the standard deviation 
is larger than a specified upper bound. This bound is a function of the wind tunnel 
facility configuration and data aquisition hardware. An upper bound of 0.01 (perturbation 
velocity) has been selected for the NASA Ames 12ft Pressure Wind Tunnel. 

CHECK 4 : WALL PRESSURE PORT NUMBER 

The total number of wall pressure ports used for the least squares fit of the wall 
signature difference is computed after CHECKS 1 TO 3 are applied. No wall interference 
corrections are computed if the total number of wall pressure ports is smaller than a 
specified lower limit. The lower limit of wall pressure ports is a function of the wind 
tunnel facility. A number of 60 wall pressure ports has been selected as the lower limit of 
the NASA Ames 12ft Pressure Wind Tunnel. 
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APPENDIX 12 


LINEAR INTERPOLATION OF WALL PRESSURE PORT CALIBRATION 

The Wall Signature Method uses a linear least squares fit of the difference between 
real-time wall signature and corresponding wall pressure port calibration to determine 
blockage effects of a test article. A value of the wall pressure port calibration has to be 
found that matches real-time test conditions. 

The wall pressure port calibration is usually a function of several independent calibra- 
tion variables. It can be a function of the total pressure pr and Mach number M at some 
wind tunnel reference station. The wall pressure port calibration can also be a function of 
the support system kinematics. In this case calibration variables describing the position 
of the support system during the calibration have to be introduced. 

An interpolation algorithm has been developed for the WICS software package that 
uses linear interpolation to determine the wall pressure port calibration as a function of 
real-time test conditions. This algorithm fulfills reliability, performance, and accuracy 
requirements of WICS. 

The interpolation algorithm is applicable to up to four independent calibration vari- 
ables. As an example, a detailed description of this interpolation algorithm for two indepen- 
dent calibration variables will be presented below. A detailed description of the algorithm 
for three and four independent variables is beyond the scope of the WICS Theory Guide. 

In general, a value of the wall pressure port calibration C(X*, Y*) has to be found for 
real-time test conditions defined by, e.g., two independent variables X* , Y* . It is assumed 
that wall pressure ports are calibrated for discrete combinations of these variables. A 
total of “n” discrete values of the first calibration variable Xi were selected. The second 
calibration variable Y{ j was changed by keeping the first calibration variable Xi constant. 
A total of discrete values of the second calibration variable Yij were selected. 

Therefore we get the following set of discrete wall pressure port calibrations : C(Xi , Yij ) 
for 1 < i < n and 1 < j < m(i) . 

The interpolation of the calibration is done in two steps. At first, calibrations are 
interpolated for the second calibration variable Y* . The first calibration variable Xi 
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is kept constant. Comparing Y* and Yij only one of three possible cases applies (see 
Fig. 41a) : 


Case 1 : Y* < 


Case 2 : Y {J < K* < Y iJ+1 


C(Xi) = C(Xi,Yij ) + 


C(X,) = C(Xi,Yi,i) 


cyww) - c(Xj,Y,j) 

J-H — 


( 12 . 1 ) 


[ y - Y,J ] (12.2) 


Case 3 : Y* > Y iiTn(i) 


C(Xi) = C(Xi,Y iMi) ) 


(12.3) 


Finally, calibrations are interpolated for the calibration variable X* by using the first 
calibration variable Xi and interpolated calibrations C(Xi) . Again, comparing X* and 
Xi , only one of three cases applies (see Fig. 41b) : 

Case 1 : X* < X x 

C(X*,y*) = C(Xx) (12.4) 

Case 2 : Xi < X* < X i+1 

C{X\Y *) = C{Xi) + — -> +1 ) ~ y (A ~ •[**-*] (12.5) 

•A* +1 - A,- L J 

Case 3 : X* > X n 

C(X*,Y*) = C(X n ) (12.6) 


Basic ideas and elements of the interpolation algorithm can easily be extended to three 
and four independent calibration variables. 

Support system wall interference corrections at a reference point grid are required if 
the Wall Interference Correction System of the 12ft PWT is applied to a fullspan model 
test configuration (see Fig. 13a) . These corrections are known as a function of calibra- 
tion variables. A similar linear interpolation algorithm may be used to interpolate the 
corrections for real-time test conditions. 
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APPENDIX 13 


SINGULAR VALUE DECOMPOSITION 

The numerical method of choice for solving linear least-squares problems, i.e. Eq. (13a) 
or Eq. (16a), is the Singular Value Decomposition (SVD) technique, [12]. This robust 
numerical technique was selected for the WICS software. The application of the SVD 
technique to Eqs. (13a), (16a) can be summarized in three steps : 

Step 1 : Rewrite Eq. (13a) or Eq. (16a) as an overdetermined linear system : 

A mX 2 ' X2xl = Bmxl (13.1) 

Step 2 : Write matrix A as the product of a column-orthogonal matrix U, a diagonal 
matrix W with positive or zero elements, and the transpose of an orthogonal matrix V . 
Matrices U, W, and V are found using the SVD algorithm. We get : 

Amx2 = U mX 2 ‘ V^2 x 2 ' v'^2x 2 (13.2) 

Step 3 : Determine the solution of the least squares problem by computing the fol- 
lowing matrix product : 

X 2 xl = V 2 x2 * W -1 2x2 * U T 2 xm ■ B mx l (13.3) 

It is important to identify elements of the diagonal matrix W, i.e. singular values of 

W, that are small. These singular values and their reciprocal will be set to zero if the 

least squares problem is ill-conditioned. 
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APPENDIX 14 


SUPPORT SYSTEM KINEMATICS 

In general, it is recommended to compute real-time coordinates of singularities that 
represent the support system or the test article if WICS is used to predict wall interference 
corrections in the 12ft PWT. This reduces the standard deviation of the least squares fit of 
the wall pressure signature and improves the accuracy of computed blockage corrections. 

The real-time coordinates of singularities are a function of the kinematics, i.e. the 
movement, of the support system or the test article. Three different types of support 
systems are presently available for testing in the 12ft PWT, i.e. the Ames Bipod, the High 
Angle of Attack Sting, and the Image Plane. Equations describing the movement of these 
support systems and of the test article in the tunnel coordinate system (see Figs. 17a, 17b) 
have to be derived. 

Ames Bipod / High Angle of Attack Sting 

The kinematics of the Ames Bipod and of the High Angle of Attack Sting in the 
tunnel coordinate system are essentially identical, if the roll angle of the High Angle of 
Attack Sting is kept at <p = 0.0° (the more complex case of a non-zero roll angle is 
discussed in Appendix 15) . The kinematics may be described by one rotation about the 
axis x* = 120.71 [ft] , z = 0.0 [ft] and one rotation about the axis x* = 120.71 [ft] , 
y = 0.0 [ft] . Angles a and (3 depicted in Fig. 42a are independent variables that 

describe the motion. For a roll angle <p = 0.0°, angle a may be approximated by the pitch 
angle of the wind tunnel model and angle {3 may be approximated by the sideslip angle 
of the wind tunnel model. Figures 42b and 42c show the connection between angle ct, 
angle /3, and singularity coordinates. It is assumed that point P\(x\,y\,z\) describes the 
initial location of a singularity. After a first rotation about the axis x* = 120.71 [/<] , 
z = 0.0 [ft] the singularity moves to point -P 2 (£ 2 , 2 / 2 > 2^)- After a second rotation about 
the axis x, = 120.71 [ft] , y = 0.0 [ft] the singularity moves to point ^ 3 (^ 3 , 2 / 3 3 ^ 3 ) - 
The final position of the singularity at point P3 is known if coordinates X3, j/3, 23 are given 
as a function of initial coordinates the pitch Jingle a, and the sideslip angle (3 . 
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Coordinates of point P$ can be derived by considering triangles depicted in Fig. 42b : 


cos 7 


X\ — x ♦ 
*1 


sin 7 


z i 


cos (7 — a) 
sin (7 — a) = 


X2 ~ x* 
R 1 


_f2_ _£3 


The following trigonometric formulas axe known : 


cos (7 — a) = cos 7 • cos a + sin 7 * sin a 

sin (7 — a) = sin 7 • cos a — cos 7 • sin a 

Combining Eqs. (14.1), (14.2), (14.3), (14.5) we get : 


X 2 = x* + (®i — X*) • cos a + z 1 • sin a 

Combining Eqs. (14.1), (14. 2), (14.4), (14. 6) we get : 

z 2 = Z 3 = zi • cos a — (n — x„) • sin a 

Considering triangles depicted in Fig. 42c we get : 

X2 — X* 


COS 7 




sin 7 


Vi 

R 2 


J/2_ 

#2 


cos (7 + /?) = 

sin (7 + /?) = 

The following trigonometric formulas are known : 


X3 — x* 
~R2~ 

J/3 

#2~ 


cos (7 + /?) = cos 7 • cos (3 — sin 7 • sin (3 

sin (7 + (3) = sin 7 • cos (3 + cos 7 • sin /? 


(14.1) 

(14.2) 

(14.3) 

(14.4) 

(14.5) 

(14.6) 

(14.7a) 

(14.76) 

(14.8) 

(14.9) 

(14.10) 

(14.11) 

(14.12) 

(14.13) 
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Combining Eqs. (14.8), (14.9), (14.10), (14.12) we get : 


i 3 = x* + (x 2 — x*) • cos P — yi • sin (3 (14.14a) 

Combining Eqs. (14. 8), (14.9), (14. 11), (14. 13) we get : 

t /3 = 3/1 • cos P + (x 2 — x*) • sin P (14.146) 

The calculation of coordinates 13 , 1 / 3 , Z 3 can be summarized as follows : Calculate X 3 
using Eqs. (14.7a) and (14.14a) . Calculate 3/3 using Eqs. (14.7a) and (14.14b) . Calculate 
Z 3 using Eq. (14.7b) . 


Image Plane 

The movement of a semispan model mounted on the Image Plane can be described by 
a single rotation about the semispan model pitch axis (x, = 120.71 [/t] , y = 0.0 [ft] ). 
The pitch angle a is the independent variable that describes the motion (see Fig. 42d). 
Figures 42e and 42f show the connection between pitch angle and singularity coordinates 
for a left and right wing semispan model. It is assumed that point Pi(x\,yi,zi) describes 
the initial location of a singularity. After a rotation about the pitch axis the singularity 
moves to the final position at point Pi{xi, 3 / 2 , zz)- The coordinate z\ of the singularity 
does not change in this case, i.e. zi = z\ . The final position of the singularity at point Pi 
is known if coordinates X 2 , 3/2 are given as a function of initial coordinates x\,y\ and the 
pitch angle a . 

Coordinates of point Pi can be derived for a left wing semispan model by considering 
triangles depicted in Fig. 42e : 


cos 7 = 

X\ — £* 

R 

(14.15) 

sin 7 = 

y 1 

R 

(14.16) 

cos (7 — a) = 

Xl - X, 

R 

(14.17) 

sin (7 — a) 

yi 

R 

(14.18) 
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Combining Eqs. (14. 5), (14. 15), (14.16), (14.17) we get : 

xz = x„ + (xi — x*)-cosa + yi • sin a (14.19a) 

Combining Eqs. (14.6), (14. 15), (14. 16), (14.18) we get : 

2/2 — 2/i • cos a — (xi — x*) • sin a (14.196) 

Figure 42f shows initial and final coordinates for a right wing semispan model. In this 
case coordinates of point Pi can be derived by simply replacing pitch angle a by —a in 
Eqs. (14. 19a), (14. 19b) and noting that cos(— a) = cos(a) and sin(— a) = — sin(a) . We 
then get : 

xi = x* + (xi — x*)-cosa — y\ • sin a (14.20a) 

2/2 = 2/1 • cos a + (®i — x*) ■ sin a (14.206) 
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APPENDIX 15 


HIGH ANGLE OF ATTACK STING KINEMATICS 


The kinematics of the High Angle of Attack Sting (HAA) for a roll angle of <p = 0.0° 
is described in detail in Appendix 14 . In Appendix 14 the reasonable assumption is made 
that the pitch angle a* of the HAA may be approximated by the angle of attack a of the 
test article and that the yaw angle /?* of the HAA may be approximated by the sideslip 
angle (3 of the test article. The line doublet orientation angle (p+ of the test article, defined 
as the angle between the positive z-direction of the test section fixed coordinate system 
(see Fig. 42a) and the lift force vector, is also equal to 0.0° for a roll angle <p = 0.0° . These 
assumptions, however, are invalid if the roll angle ip of the HAA is not equal to 0.0° . 

Equations are derived in this Appendix that use the angle of attack a, sideslip angle (3 , 
and roll angle <p of the test article to (i) compute the pitch angle a* and the yaw angle j3* 
of the HAA, to (ii) compute the location of test article singularities and reference points, 
and to (iii) compute the orientation angle of test article line doublets. 

(i) PITCH ANGLE a* AND YAW ANGLE /?* OF THE HAA : The pitch and yaw 
angle of the HAA may be computed in several steps assuming that angle of attack a, 
sideslip angle /?, and roll angle <p of the test article are known. In a first step, the unit 
wind vector is expressed in the model coordinate system defined by unit vectors i , 
j , and k . Using Fig. 43a it can easily be seen that 

( cos a • cos (3 \ 

- sin (3 I (15.1) 

sin a • cos (3 J 


In the second step, it is necessary to reverse the roll angle rotation. Assuming that 
the roll axis of the test article is identical with unit vector i (see Fig. 43b), we get for 
the rotated unit vectors / , J , and K : 



J = 



(15.2a) 


(15.26) 
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It = 


0 

— sin <p 
cos f 


(15.2c) 


Now it is possible to compute the yaw angle /?* of the HAA by using a scalar product 
and a right-angled triangle depicted in Fig. 43c . The yaw angle /?* is defined as the angle 
between unit wind vector W and the plane spanned by unit vectors ~f and it . Using 
the projection [Woo o J ] J of the unit wind vector Woo in the direction of unit vector 
J and the sign convention that /?* is positive if the unit wind vector is coming from the 
right wing, we get : 

sin /?* = - Woo o J (15.3a) 


or 


/?* = arcsin 


~ Woo o J 


_ JL < < JL 

2 “ y ~ 2 


(15.36) 


where the scalar product W Q 0 o J is computed by using Eqs. (15.1) and (15.2b) : 


Woo 0 J — — sin (3 • cos (f + sin a • cos f3 • sin 


(15.3c) 


Similarly, the pitch angle a* of the HAA may be computed by using a scalar product 
and a right-angled triangle depicted in Fig. 43c . The pitch angle a* is defined as the angle 


between unit vector I and the projection of the wind vector Woo on the plane spanned 
by unit vectors / and K . Using the projection of vector Woo — [ Woo o J ] ~J in the 
direction of unit vector it we get : 


sm a* = 


( W TO - [ Woo o J ] J ) o K 


Woo O J 


(15.4a) 


Knowing that 7 1 it, i.e. (Woo - [w2o7]7)olt = V/Z olt , we finally get 


a* = arcsin 


Woo o K 


1 - 


Woo o J 


T 

~2 


7T 


- - T 


(15.46) 


Scalar products in Eq. (15.4b) Eire computed by using Eqs. (15.1), (15.2b), and (15.2c) . 
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(ii) SINGULARITY AND REFERENCE POINT COORDINATES : It is now possible 
to compute the location of singularities and reference points of a test article that is mounted 
on the HAA by using Eqs. (15.3b), (15.4b) and Eqs. (14.7a), (14.7b), (14.14a), and (14.14b) 
of Appendix 14 . This is done as follows : 

( 1 .) Apply roll angle <p to the initial singularity or reference point location given by point 
Po(x<>,yo,Zo) • Using Fig. 43d it is possible to compute the new location of the singularity 
or reference point at point Pi(®i , j/i , zi) as : 


x\ = x 0 (15.5a) 

j/i = R 0 ■ cos [ fi - <p ] (15.56) 

Z\ = R 0 • sin [ ft — tp ] (15.5c) 

Combining Eqs. (15.5b), (15.5c), (14.5), (14.6) and knowning that cos p. = y 0 /Ro and 
sin fi = Zo/Ro we finally get : 

y\ = 2/o • cos <p + z 0 - sin ip (15.5d) 

zi = z 0 • cos <p — y 0 ■ sin <p (15.5e) 


(2.) Rewrite Eqs. (14.7a), (14.7b), (14.14a), (14.14b) by replacing a, f3 with a*, /?, : 


X 2 = X, 

+ (xi — x*) ■ cos a* + z\ • sin a * 

(15.6a) 

X3 = X* 

+ (x 2 - x*) • cos /?* — y\ • sin /?, 

(15.66) 

yz = 

2/i • cos (3+ + (x 2 - x,) • sin (3* 

(15.6c) 

Z 3 = 

zi • cos a» — (xi — x*) • sin a* 

(15. 6d) 

The calculation of coordinates X3, 2/3, z$ may be summarized as follows : Calculate X3 
using Eqs. (15.6a) and (15.6b) . Calculate 2/3 using Eqs. (15.6a) and (15.6c) . Calculate 


Z 3 using Eq. (15.6d) . 

(iii) LINE DOUBLET ORIENTATION ANGLE <p . : In general, the fine doublet 
orientation angle <p* of the test article (identical with angle r in Appendix 4) is defined as 
the angle between the positive z-direction of the test section fixed coordinate system (see 
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Fig. 42a) and the lift force vector. Angle p* can easily be computed using a unit vector 
that is parallel to the positive z-direction of the test section coordinate system and a unit 
vector that is parallel to the lift force vector. 

A unit vector Z that is parallel to the positive z-direction of the test section co- 
ordinate system may be obtained by rotating unit vector K in the I - K plane (see 
Figs. 43c, 43e) . The selected rotation angle is equal to the pitch angle a* of the HAA 
given in Eq. (15.4b) . Using Fig. 43e we then get : 


Z = 


K — tan a* • I 
K — tan a* - / 


(15.7a) 


where 


it - 


tan a* 


/ — tan a* \ 


/ = I — sin ip 1 

(15.76) 

\ COS ip J 


1 j = \/ 1 -J- tan 2 a* 

(15.7c) 


A unit vector W that is parallel to the lift force vector may be obtained by computing 
the cross-product Woo x j using Eqs. (15.1) and knowing that ~j = ( 0 ; 1 ; 0 ) (see 
also Fig. 43a) . We then get : 


W = 


Woo X j 


Woo X j 



(15.8) 


Finally, the line doublet orientation angle <p* may be obtained by computing the scalar 
product of unit vectors W and Z . We get : 


vt o Z = cos <p+ 

Combining Eqs. (15.7a), (15.7b), (15.7c), (15.8), and (15.9) we obtain : 


<P * = 


arccos o Z ; <p < tt 


<p * = 2 i r — arccos ^ W o Z 


where 


^ ; V > tt 

-=* sin a • tan a* + cos a • cos <p 

W o Z = 




(15.9) 


(15.10a) 

(15.106) 

(15.10c) 


86 



APPENDIX 16 


PITCHING MOMENT LEVER ARM 


In Appendix 5 equations axe derived that compute the singularity strength of a line 
doublet of a wind tunnel model by using the measured lift force and pitching moment. 
The pitching moment lever arm is approximated in Eq. (5.2) as the difference between 
the x-coordinate of the pitching moment reference axis and the x-coordinate of the line 
doublet starting point (see Fig.38) . This approximation, however, is not valid anymore if 
a wind tunnel model is tested at a high angle of attack. More general equations are derived 
in this Appendix that use the pitching moment axis and the lift force axis to compute the 
pitching moment lever arm of a line doublet. 

In general, the pitching moment lever arm of a line doublet is defined as the distance 
between the pitching moment axis and the lift force axis of a line doublet (see Fig. 44) . 
Therefore a mathematical description of each axis has to be found in order to compute the 
pitching moment lever arm. 

(i) PITCHING MOMENT AXIS EQUATION : The pitching moment axis equation 
may be found by using the initial coordinates of the balance center in the tunnel coordinate 
system (see Fig. 17a, 17b). Assuming that x mr , y mr ,Zmr are the initial coordinates of the 
balance center (i.e. angle of attack, sideslip angle, and roll angle of the test article are 
equal to zero), it is possible to introduce an initial unit pitching moment axis vector A as 
follows : 


A 



(16.1 a) 


where 

( Xmr \ 

A 1 = I J/mr I (16.16) 

\ z mr ) 

For a fullspan model point vector A 2 is defined as : 

Xrnr \ 

ymr+l.Oj (16.1C) 

z mr / 
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For a left wing semispan model point vector A 2 is defined as : 

A2 = 



(16.1 d) 


For a right wing semispan model point vector A 2 is defined as 

A2 = 



(16. le) 

In the next step, equations describing the kinematics, i.e. the movement of the model 
support system, have to be applied to the points described by vector A\ and A 2 (see 
Appendix 14 and 15). Point vector moves to point vector (if and point vector A 2 
moves to point vector C* . After points A 1 , A 2 are moved to their new location, it is 
possible to describe the unit pitching moment axis vector D as : 


d = c* 


(16.2) 


A point P on the pitching moment axis may now be described by using vectors C?, ~D , 
and an independent variable A : 


P( A) = C + A - D 


(16.3) 


(ii) LIFT FORCE AXIS EQUATION : The lift force axis equation of a line doublet 
may be found by using the tunnel coordinates of the line doublet starting point in com- 
bination with the line doublet orientation angle (see also Appendix 4 and 15) . After 
the application of equations describing the support system kinematics, the line doublet 
starting point has moved to a location described by point vector E . The line doublet 
orientation is given by a unit lift vector F in the tunnel coordinate system : 

0 ) 

sin tp* I (16.4) 

cos <p* J 

For a left wing semispan model p>+ = 7r/2, for a right wing semispan model tp „ = — tt/2 . 
Similar to Eq. (16.3), a point Q on the lift force axis may be descibed by using vectors 
E , F , and an independent variable p : 


Q(p) = E + p- F 


(16.5) 
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(iii) PITCHING MOMENT LEVER ARM : The pitching moment lever arm of the 
selected line doublet is found by minimizing the square of the distance of points located on 
each axis (see Fig. 44). The square of the distance of two selected points may be computed 
using a scalar product : 

D 2 (\n) = [ f?(A) - ] o [ f?(A) - 3o<) ] (16.6) 


Using Eqs. (16.3), (16.5) and applying the product rule to Eq. (16.6) we get : 

dD\ A,/i) 


d A 

d £> 2 ( A,ji) 
d n 


= 2 • [ P (A) - Q(p)) o D 


= 2 • [ P (A) - <}(/i)] o [- F } 


(16.7a) 


(16.76) 


The variables A mtn and of the two closest points are found by setting Eqs. (16.7a), 
(16.7b) to zero. We get : 


[ P (Amin ) - Q(/lmin)] O D = 0 
[ P(A m ,n) - Q(/W>) ] O F = 0 


(16.8a) 


(16.86) 


Using Eqs. (16.3), (16.5), it is possible to write Eqs. (16.8a), (16.8b) as a 2 x 2 system of 
linear equations : 

[~D o ~D ] • A min + [ - D o ~F ] • fJLmin = [ ~E - ] o ~D (16.9a) 

[D O P ] . X min + [ - P O P = [ E - ~6 } O P (16.96) 

In general, D L F or D o F = 0. Therefore, we get for A mtn and n m in of the two 

closest points : 

. [ P - 15 } o D 


Umin — 


D O D 

[ ~E - ] o ~F 


(16.10a) 

(16.106) 


- F o F 

Finally, using Eq. (16.6), we get for the pitching moment lever arm : 


7^(A m ^ n } flrnin ) — \*[ P ( A m j n ) Q (/^mtn ) ] ° [ P(Amtn) Q (^rntn) ] (16.11) 
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The following two conditions may be used to determine the correct sign of the pitching 
moment caused by a line doublet : 

[ D X y ] O [ P(X min ) - Q(/W) ] > 0 =* “+" 

[ D x F ] o [ P(A min ) - Q(^ in ) ] < 0 

These conditions ensure that the pitching moment is positive if the lift force axis of a 
line doublet is upstream of the pitching moment axis (see also Eq.(5.2) and Fig. 38) . 
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APPENDIX 17 


INCLINATION OF FORCE AND MOMENT VECTORS 


Corrections to the lift, drag, rolling moment, and yawing moment coefficient are caused 
by the wall interference induced inclination of corresponding force and moment vectors. 
These four corrections may be estimated by describing the inclination of the force and 
moment vectors using the induced mean angle of attack correction a; at the 1/4-chord 
line of the wing (see also Ref. [18], p.17). 

LIFT AND DRAG COEFFICIENT CORRECTION 

The lift and drag coefficient correction caused by the inclination of the lift and drag 
force may be estimated by using the resultant of the uncorrected lift and drag force vectors. 
Using Fig. 45a we get for the resultant vector in the x-z coordinate system : 

S’ = 0 + V = (^) (17.1) 


where D' is the magnitude of the uncorrected drag force and V is the magnitude of the 
uncorrected lift force experienced by the wind tunnel model. 

Unit vectors i , j in the free-stream direction and perpendicular to the free-stream 
direction of the wind tunnel flow field have to be specified as they are used to define the 
lift and drag force. These unit vectors axe a function of the mean angle of attack correction 
ai at the 1 /4-chord of the wing. They may be written in the x-z coordinate system as : 


J 



(17.2a) 

(17.26) 


The mean angle of attack correction a* may be computed by calculating the local 
angle of attack correction at equally spaced reference points along the 1 /4-chord line of 
the wing. We then get : 


1/4 — chord line 


Oti — 


Ef=i 


Ef= i Mi) 


(17.2c) 
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where N is the number of reference points, w is the weighting factor of a reference point (a 
function of the lift distribution of the wing), and a, is the local angle of attack correction 


at a reference point. Now it is possible to compute the lift and drag force vectors 

as : 




• , • 





L = L j 


R o j 

j 

(17.3a) 




r 




D = D i 


R o i 

i 

(17.36) 

where 






L 

— R o ~j = 

- D' • 

sin Si + 

T f ^ 

L • cos ot{ 

(17.4a) 

D 

= R o~i = 

D' • COS CCi + 

ft • 

L • sin (*i 

(17.46) 


The lift and drag coefficient correction may now be written as : 


A CL = 

CL - c' L = 

L - V 

qoo ■ s 

(17.5a) 

A C£> = 

CD — c' D = 

D-D' 
?oo * S 

(17.56) 


where is the free-stream dynamic pressure and S is the reference area of the wind 
tunnel model. A first order approximation of the free-stream dynamic pressure may be 
obtained by using Eq. (1.15) in Appendix 1. Using the test section reference density, 
velocity, and Mach number we get : 


9oo — ?re/ • [ 1 + ( 2 — ) • c ] 

= -?f- ■ U?', • [1 + (2 - Ml, ) • ?] 

The mean blockage factor e has to be computed as a. weighted average of local blockage 
factors e that are computed at a set of user selected reference points. We get : 

_ _ SjLi w (j) •*(;') 

£f=i t v(j) 

where N is the number of reference points, w is the weighting factor of a reference point, 
and e is the local blockage factor at a reference point. 




92 



Important remark: The reference points used to compute the mean blockage factor 
in Eq. (17.6b) are not necessarily identical with the reference points used to compute the 
mean angle of attack correction in Eq. (17.2c) ! 

Combining Eqs. (17.4a), (17.4b), (17.5a), and (17.5b) we get for the lift and drag 
coefficient correction : 


A cl 


A CD 


— D' • sin cti + V • [ cos Si — 1 ] 
9oo ■ S 

D r • [ cos ati — 1 ] + V ■ sin Si 
9c© ' S 


(17.7 a) 
(17.76) 


Finally, corrections have to be added to the uncorrected lift and drag coefficients. 
Then we get : 

cl,c = CL,unc -I - A cl (17.8a) 

CD,c ~ CD, tine + A c D (17.86) 


ROLLING AND YAWING MOMENT COEFFICIENT CORRECTION 

The rolling and yawing moment coefficient correction caused by the inclination of 
the rolling and yawing moment vector may be computed by using the resultant of the 
uncorrected rolling and yawing moment vectors. Using Fig. 45b we get for the resultant 
vector in the x-z coordinate system : 

f = # + Y' = (~y,) (17.9) 

where R' is the magnitude of the uncorrected rolling moment and Y' is the magnitude of 
the uncorrected yawing moment experienced by the wind t unn el model. 

Similar to Eqs. (17.3a) and (17.3b), it is possible to combine Eqs. (17.2a), (17. 2b) with 
Eq. (17.9) . Then, we get for the rolling and yawing moment vectors : 


where 


R = — R ~i 

= 7o7 

i 

(17.10a) 

Y = - Y 7 

= 7 o 7 

. 

3 

(17.106) 

= — 7 o i = 

R' • cos Si + 

Y 1 • sin Si 

(17.11a) 
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Y = — if o j = — R' • sin cfi + Y’ - cos Si (17.116) 

The rolling moment and yawing moment coefficient correction may now be written as : 


A ci 


R - R' 

goo • S • 6/2 


(17.12a) 


Ac n — c n c n 


Y_-_Y^_ 

qoo-S- 6/2 


(17.126) 


where goo is the free-stream dynamic pressure, S is the reference area of the wind tunnel 
model, and 6 is the wing span. 

Combining Eqs. (17.11a), (17.11b), (17.12a), and (17.12b) we get for the rolling and 
yawing moment coefficient correction : 


A cj 


R' ■ [ cos Si — 1 ] + Y' ■ sin Si 

q<x • S • 6/2 


(17.13a) 


A _ _ — R' ■ sin Si + Y' ■ [ cos <5,- — 1 ] 


(17.136) 


?oo • S • 6/2 

Finally, corrections have to be added to the uncorrected rolling and yawing moment 
coefficients. Then we get : 


C l,c — Cl,unc “f* Acj 

(17.14a) 

C n,c — Cn,unc "h Ac n 

(17.146) 
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APPENDIX 18 


PITCHING MOMENT COEFFICIENT CORRECTION 


Wall interference induced angle of attack corrections are not constant in the test 
section of a wind tunnel (see Fig. 24c). Therefore, it is necessary to compute a pitching 
moment coefficient correction that will account for the difference between the mean and 
local angle of attack correction in the spanwise and chordwise direction of the wing of the 
test article. 

In general, the corrected pitching moment coefficient cm,c of the complete test article 
configuration may be written as : 

c Af ,c ~ C Af t tm c + Acmi + Acm 2 ( 18 . 1 ) 


where CM.unc is the uncorrected pitching moment coefficient, Acm i is the pitching moment 
coefficient correction due to the difference between the mean and local angle of attack cor- 
rection in the spanwise direction of the wing, and Acm 2 is the pitching moment coefficient 
correction due to wall interference induced streamline curvature in the chordwise direction 
of the wing. 

PITCHING MOMENT CORRECTION Acmi : 


A pitching moment coefficient correction Acmi is due to the difference between the 
mean and local angle of attack correction in the spanwise direction of the wing. It may 
also be interpreted as a correction resulting from the lateral shift of the spanwise center of 
lift caused by wall interference effects. 

The pitching moment coefficient correction Acmi may be computed by using Eq. (3.58) 
given in Ref. [10], which is essentially identical with Eq. (34) in Ref. [18] . Using the no- 
tation introduced in the WICS Theory Guide we get : 

Acmi = F 'f l^o~ ’ d ( 18 - 2a ) 


180 
F = 


(7t/2) • A 2 • (dcL/da) • tan ^> 0.25 


7T-A + 2 -(dc L /da) (18,26) 

where A is the aspect ratio of the wing, dc^/da is the lift curve slope in [ 1 /rad], 920.25 is 
the sweep angle of the 1/4-chord line of the wing, 0 ^( 0 . 75, y) is the local angle of attack 


95 



correction in [deg] along the 3/4-chord line of the wing, c*7 is the angle of attack correction 
in [deg] of the wing, w w is a weighting factor derived from the lift distribution (see also 
Eq. (19.2c) in Appendix 19), y is the spanwise coordinate, and 6 is the span of the wing. 

In general, Simpson’s Rule may be applied to Eq. (18.2a) if Acm i has to be determined 
for a given wing. However, it is also possible to simply discretize Eq. (18.2a) as long as the 
local angle of attack correction a t -(0.75, y) is computed at a large number of equally spaced 
points along the 3/4-chord line of the wing. Assuming that a total number of N equally 
spaced reference points were selected along the semispan of the wing (i.e. 6/2 = N ■ Ay), 
we get the approximation : 

f 2 y \ 2 1 

d {— ) K T Ay - i r ( 18 - 3 -) 

Using the abbreviation w w (y(j)) = w w (j), Eq. (18.2a) may be discretized as : 

N 

Acm\ « F ’ Yl "j^o" ' a *'(°- 75 ’-?) - “i] * W w(j) * — (18.36) 

j= 1 

After rearranging terms we get : 

F • n 1 N 

Acmi ~ ~b~ 90 ' ~W ' [ a <(°- 7 M “ “i] ■ W M • y(j) (18.3c) 

i=i 

Sometimes it is necessary to distribute a total number of M equally spaced reference 
points along the complete span of the wing (i.e. 6 = M ■ Ay). Then, using the absolute 
value of y(j), Eq. (18.3c) becomes : 

Acmi ~ 1 • 90 ’ m ' ^ [ a *(°- 75 >i) - a *-J • W *U) • ly(i)l (18.4) 

i=i 

WICS uses Eq. (18.4) to determine the pitching moment coefficient correction Acjvfi 
that is caused by the wall interference induced shift of the center of pressure of the wing. 
PITCHING MOMENT CORRECTION Ac M 2 : 

A pitching moment coefficient correction Ac\f% due to the streamline curvature along 
the chord of the wing has to be computed. Classical thtory has shown that the magnitude 
of this correction may be approximated by the camber of a two-dimensional circularly 
cambered airfoil (see Ref. [18], p.17/18 and Ref. [28], p.91) . 
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In general, the angle of attack correction a,- is not constant along the chord of a wing. 
This phenomena may be described as a wall interference induced streamline curvature that 
is experienced by the wing. A constant streamline curvature may be found if we assume 
a linear change of the angle of attack correction along the chord (see Fig. 46a) . Then, a 
linear function 6a t - may be defined as (see Fig. 46b) : 

= [ a i(0 ~ a* (0.50) ] • ^ [rad] (18.5) 

The linear function 6a,- may be related to a two-dimensional circularly cambered 
airfoil, [28], because its derivative with respect to £ is the streamline curvature. The pitching 
moment coefficient about the 1/4-chord point of a circularly cambered airfoil is given in 
Ref. [28] as : 

cjvf(0.25) = - 7T • — (18.6) 

c 

where f/c is the camber of the airfoil (see Fig. 47). The camber is related to the slope of 
the circular arc at the 1/4 and 3/4-chord point. The slope is related to the upwash along 
the chord. Then, using Fig. 47, we get : 

„ „ 1 
tan p « p = — • 

2 

Using triangle ABC in Fig. 47 we also know that 

tan ft « j3 = — — (18.76) 

c 

Combining Eqs. (18.7a), (18.7b) we get for the camber as a function of the streamline 
curvature : 

— = 4" * [ 6a, (0.25) - 6a*(0.75) ] (18.8) 

c 4 

Using Eq. (18.5) we also know : 

6a,(0.25) = [ a,(0.25) - a,-(0.50) ] • -^- [rad] (18.9a) 

180 

<5ai(0.75) = [ a<(0.75) - a,-(0.50) ] • ^ [rad] (18.96) 


d z 

IT 


d z 

~dj 


0.25 \ ' / 0.75 

6a, (0.25) - 6a,-(0.75) ] 


(18.7a) 
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Combining Eqs. (18.6), (18.8), (18.9a), (18.9b) we finally get for the pitching moment 
coefficient of the two-dimensional circularly cambered airfoil as a function of the difference 
between the angle of attack correction on the 1/4 and 3/4-chord point : 


cm(0.25) = [ <*(0.75) - <*(0.25) ] • ~ (18.10) 

In the next step, it is necessary to integrate the pitching moment coefficient in the 
spanwise direction of the wing. Assuming that the sweep angle <p of the wing is 0.0° (see 
Fig. 48), the local pitching moment coefficient of a wing segment of width dy about the 
1/4-chord line may be written as 


CAf(0.25,y) 


d Pjy) _ 

goo ■ c(y) • dS(y) 


(18.11) 


where d P(y) is the infinitesimal pitching moment acting on the wing segment, <joo is the 
free-stream dynamic pressure seen by the wing segment, c(y) is the local chord, and d S(y ) 
is the infinitesimal area of the wing segment. The infinitesimal area d S(y) of the wing 
segment is a trapezoid and may therefore be expressed as : 


d S(y) = c(y) • d y 


(18.12) 


Introducing the mean geometric chord as a function of the wing reference area 5 and the 
wing span b we know : 

S , 

c = — (18.13) 

In the next step we combine Eqs. (18.11), (18.12), divide both sides of the resulting equation 
by the mean geometric chord c and the mean aerodynamic chord d, introduce Eq. (18.13) 
on the right hand side of the equation, and rearrange terms. Then, we get : 

c 2 (y) 1 __ d P(y ) 


caf(0.25,y) 


dy = 


(18.14) 


c'-c b " qoo . d • S 

The pitching moment coefficient correction Acm 2 for 0.0° sweep angle may be obtained 


by integrating both sides of Eq. (18.14) over the whole wing span : 


r b / 2 

Acm2 = / CA/(0.25,y) 

J-b/ 2 


c 2 (y) d y 
d • c b 


■ i: 


dP(y) 


Qoo * d * S 


(18.15) 


98 



The pitching moment coefficient of a swept wing may be approximated (see Ref. [29]) 
by using the product of the pitching moment coefficient of the unswept wing and the 
cosine of an average sweep angle <p. For a circular cambered airfoil the sweep angle at the 
1 /2-chord is selected. We get : 


AcAf 2 (9 sl 7^ 0.0°) = cos ipo.so • Ac M 2 (<P = 0.0°) 


(18.16) 


Finally, combining Eqs. (18.10), (18.15), (18.16) we get for the pitching moment coefficient 
correction Acm2 due to streamline curvature : 

7T • COS 9^0.50 f b/2 r , n -- ^ /n or; M x ^(.V) d V 

[ a t (0.75, y) - a f (0.25,y) ] • : 

- 6/2 


ACM2 = 


£>-> ~ nr (1817) 


The variable transformation /i = 2 y/6 is introduced in an effort to compare Eq. (18.17) 
with an equation reported in the literature. Assuming that the lift distribution is symmet- 
ric, Eq. (18.17) may then be written as : 


Acjif 2 = 


7T • COS ¥>0.50 


f ^ TT C^( itl 

•/ [ ai(0.75,ft) - a, -(0.25, /i) ] • — • d y, (18.18) 


Equation (18.18) agrees with Eq. (35) in Ref. [18] . 

Assuming that local Jingle of attack corrections are computed on a total of N equally 
spaced reference points along the 1 / 4-chord line of the wing and N equally spaced reference 
points along the 3/4-chord line of the wing (i.e. 6/2 = N • Ay), we get the approximation : 

2 Ay 1 


d fi 


An = 


(18.19) 


b N 

Combining Eqs. (18.18), (18.19) we get an approximation of the pitching moment coeffi- 
cient correction due to streamline curvature : 

N - ,2( 


ACM2 


7T • COS <f 0.50 


*53 _ ttf(0.25, j) ] • ^ 

i=i 


c 2 U) 


N 


(18.20a) 


Sometimes it is required to distribute a total number of M equally spaced reference 
points along the 1/4- and 3 /4-chord line of the complete span of the wing (i.e. 6 = M • Ay) . 
Then Eq. (18.20a) becomes : 

7T- COS ¥>o.50 ^ T ^ .-)! * <?<J) 


ACM2 


««(0.75,j) - «i(0.25,j) • ^ 


3 = 1 


180 d ■ c 


M 


(18.206) 
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Finally, introducing compressibility effects in Eq. (18.20b) (see Ref. [10], Eq. (3.59)), 
we get : 


A CM2 



M 

£ [ <*(0.75, i) - <*(0.25,;) ] ■ c\j) 


(18.21a) 


where 


G 


7T • COS y>0,50 w 1 

4 ■ y/T—~M^~cos^ipo^o 180 d • c 


(18.216) 


Equation (18.21a) and (18.21b) are used by WICS to determine a pitching moment 
correction due to streamline curvature. 


CALCULATION OF PITCHING MOMENT COEFFICIENT CORRECTIONS : 

WICS uses Eq. (18.4) and Eqs. (18. 21a, b) to compute pitching moment coefficient 
corrections due to the difference between the mean and local wall interference corrections 
in the spanwise and chordwise direction of the wing of a wind tunnel model. Input variables 
required for the application of these equations axe classified as follows : (1) variables that 
depend on the geometry of the test article (A, c, c', c, <p Q . 25 , po.so), (2) variables that 
depend on the lift of the wing (w w (j), dcL/da ), and (3) variables that are a result 
of the calculation of mean and local wall interference corrections (a7, a t -(j); a,- (0.25, j), 
a,-(0.75, j)). For best results it is recommended to determine the lift curve slope dcLjda 
by using the experimental cl — a curve. 
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APPENDIX 19 


ROLLING MOMENT COEFFICIENT CORRECTION 

An additional rolling moment coefficient correction Ac/ has to be found whenever 
wall interference effects are the result of a non-symmetric lift distribution on the wing. A 
non-symmetric lift distribution may be caused by deflected ailerons or non-zero sideslip 
angle of the wind tunnel model. 

In general, the corrected rolling moment coefficient c/ iC of the complete test article 
configuration may be written as : 

Cj,e = Cl,unc + Ac/ (19.1) 

where c/ itinc is the uncorrected rolling moment coefficient and Ac/ is the rolling moment 
coefficient correction due to wall interference effects. 

The calculation of the rolling moment coefficient correction is done in three steps. At 
first, it is necessary to find a reasonable non-symmetric lift distribution and fine doublet 
representation of the wing that satisfies the measured lift force, pitching moment, and 
rolling moment of the test article. Then, angle of attack corrections along the 3/4-chord 
of the wing are computed using the non-symmetric line doublet representation of the wing 
and the wall signature method. Finally, the rolling moment coefficient correction Ac/ is 
determined by integrating the angle of attack correction distribution. 

LIFT DISTRIBUTION AND LINE DOUBLET REPRESENTATION : 

The calculation of a rolling moment coefficient correction is only possible if a non- 
symmetric lift distribution of the wing is found that matches the measured lift force, 
pitching moment, and rolling moment of the test article. This lift distribution may be 
obtained in two steps by using results presented in Appendix 5. 

At first, a symmetric lift distribution of the wing is assumed. The corresponding 
circulation T* of the wing is written as (see Appendix 5 for more detail) : 

r *{y) = T w • w w (y) (19.2a) 

where T,,, is only a function of the lift force and pitching moment measurement on the test 
article and w w (y) is the weighting factor of the selected symmetric lift distribution. The 
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weighting factor w w (y) has to fulfill the condition : 



(19.26) 


The weighting factor w w (y) of an elliptic lift distribution (6 = wing span) is given as 


w, 


>(y) = 


A_ 

IT 


rrr 

2-y 1 

1 

b 


1 2 


(19.2c) 


The total hft caused by the symmetric lift distribution may be computed by integrating 
the circulation along the wing span. Using Eqs. (19.2a), (19.2b) we get : 


/• ft /2 

^(r ) = Poo ■ ■ / r (y) d y = /?oo • Uqo • 

J-b/2 


(19.3a) 


Similarly, we get for the total pitching moment caused by the symmetric lift distribu- 


tion : 


rft/2 

- ft /2 

= Poo ' Woo ' Tai ‘ 


= Poo ' «oo ' f A x(y) ■ r*(y) d y 

J-b/2 

yft/2 

/ Ax(y) • w w (y) d y 

J-b/2 


(19.36) 


where A x(y) is the pitching moment arm. 

In a second step, a weighting factor Aw is applied to the symmetric lift distribution 
such that the total rolling moment caused by the resulting non-symmetric lift distribution 
equals the measured rolling moment R. The weighting factor Aw must be introduced such 
that the total lift force and pitching moment of the non— symmetric hft distribution equals 
the lift force and pitching moment of the symmetric lift distribution. This weighting factor 
Aw may be defined as follows (see Figure 49) : 


r ~(v) 

= r *(y) • [ 1 + Aw ] ; 

A 

o 

(19.4a) 

r **(») 

= r *(y) ; 

y = 0 

(19.46) 

r *•(») 

= r*(y) • [ 1 - Aw ] ■ 

y > 0 

(19.4c) 


where T** is the circulation of the non-symmetric lift distribution of the wing. 
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Combining Eqs. (19.2a), (19.4a), (19.4b), and (19.4c) it is possible to express the 
non-symmetric lift distribution as a function of the symmetric lift distribution. Similar to 
Eq. (19.2a), we get for the circulation of the non-symmetric lift distribution : 


where 



r**(y) = T w • w w (y) 


(19.5a) 

w^(y) 

= w w (y ) • [ 1 + A w ] ; 

y < 0 

(19.56) 

w^(y) 

= Wu ,(y) ; 

y = 0 

(19.5c) 

wZ(y) 

= W w (y) • [ 1 - Aw ] ; 

y > 0 

(19. 5d) 


In the next step it is necessary to show that the lift and pitching moment caused by the 
non-symmetric lift distribution T** is identical with values obtained from the symmetric 
lift distribution T*. Integrating the circulation T** of the non-symmetric lift distribution 
along the wing span and introducing Eqs. (19.5a) we get : 

r &/2 


/ o// 

r "(y) d y 

■b/2 

/• 6/2 

= Poo * ^oo I w w (y) d y 

J- b/2 


(19.6a) 


Using Eqs. (19. 2b), (19. 5b), (19. 5c), (19. 5d) and knowing that w w (y) is symmetric to y = 0 
we get : 


rW 2 

w^{y) d y 



J-b/ 2 




r° 


r b / 2 


= / 

w^{y) dy + 

/ w^(y) d y 


J —b/2 


Jo 


f 0 


[b/2 


= 

w w (y) * [ 1 + 

Aw ] d y + / w w (y) • 

[ 1 — Aw ] d y 

J-b/2 


Jo 


fW 2 


[ yo 

[b/2 

= 

w w (y) d y + 

Aw • / w w (y) dy - 

/ w w (y) d y 

J-b/2 


\J-b/2 

/ 0 


(19.66) 


= 6 

Thus, combining Eqs. (19.6a), (19.6b) and comparing the result with Eq. (19.3a) we con- 
clude that : 

L(r*) = Poo • «oo • r„ • 6 = L(F) (19.6c) 
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Similarly, we get for the pitching moment of the non-symmetric lift distribution : 


P( r **) = Poo-Uoo- f A*(y) ■ r*(y) d y 

J-b / 2 

rW 2 

= Poo • u<x> -T w ■ I Ax(y) • w w (y) d y 

J-b/2 


(19.7a) 


Using Eqs. (19.5b), (19.5c), (19.5d) and knowing that A x(y) • w w (y) is symmetric to y = 0 
we get : 


r b/2 

Ax(y) 

-6/2 


■0 

A x(y) 

-6/2 


0 

A x(y) 

-6/2 

/ 6/2 

+ 

l 

6/2 

A x(y) 

-6/2 


+ 

Au; • 

6/2 

Ax(y) 

•6/2 



f° fb / 2 

/ Ax(y) ■ w w (y) dy + / A x(y) • w^(y) d y 

J-b/2 JO 

I 

A x(y) • w w (y) • [ 1 + Au; ] d y 

b/2 

fb/2 

+ / Az(y) • w w (y) • [ 1 — Aw ] d y 

Jo 

/2 

Ax(y) • w w {y) d y 

S/2 

+ Au; • /:, A x(y) • w w (y) dy - A x{y) 


(19.76) 


w w (y) d y 


Thus, combining Eqs. (19.7a), (19.7b) and comparing the result with Eq. (19.3b) we con- 
clude that : 


f b / 2 

P{ r**) = Poo * «oo * r«, • / A x(y) • W w iy)dy = P(r*) (19.7c) 

J-b/2 

In Eqs. (19.6c) and (19.7c) it was shown that the introduction of the weighting factor 
Atu does not change the lift force and pitching moment.. It remains to determine Au; as a 
function of the rolling moment measurement. This may be done as follows : 

The calculation of the angle of attack correction based on the w all signature method 
requires that the non-symmetric lift distribution of the wing has to be translated to a 
corresponding non-symmetric line doublet representation of the wing. We assume that a 


104 



total number of n w equally spaced line doublets axe selected along the 1/4-chord line of 
the wing. Then, the strength u^*(j/(j)) of a line doublet at the spanwise location y(i) is 
defined as (see also Eq. (5.4a) in Appendix 5) : 


°C(y(*)) = As* * r„ • w w (y(i)) 


(19.8a) 


where w w (y(i)) is given by Eqs. (19.5b), (19.5c), (19.5d) and 

6 


A s w = 


n v 


y{i) = 


-b 


+ A 


i — 


(19.86) 


(19.8c) 


The total rolling moment caused by the non-symmetric line doublet representation of 
the wing may now be written as : 


~ R = 2 “M*) ‘ »(*') 


(19.9a) 


t = l 


where L w (i) is the lift force caused by a single line doublet at the spanwise station y(i) . 
Equation (19.9a) may also be written as the sum of the contributions from the left and 
right wing. For an even number n w of line doublets we then get : 


W /2 fl *£( 

-r = 2 Lw ^ ■ + ^*(0 ■ y (*) 

* = 1 X w j 2 -f* 1 


(19.96) 


Applying the Kutta/Joukowski formula to the wing span increment As^,, introducing 
Eq. (19.5a), and using the abbreviations r**(y(i)) = r**(i) and w^(y(i)) = w w (i) we get : 


L w (i) — Poo * Uqq • As w • r ( 2 ) — Poo * a oo • A ■ 10 ^( 2 ) 


Combining Eqs. (19.5b), (19.5c), (19. 5d), (19.9b) and (19.10) we get : 


(19.10) 


-R 


n w /2 


Poo * ^00 * * r* u 


= 1 + A w ] • ^2 W w(i) • y(i ) 


t=l 


(19.11) 


+ [ 1 — Aw ] • E *M*) • y(0 

/2+1 
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The weighting factor w w (i) of the symmetric lift distribution is symmetric to the root of 
the wing (y = 0). Therefore, we know : 


n«, n w /2 

53 u>w(i) • y(i) = - 53 “’■(*) ' KO (19.12) 

i=n w /2+1 „ t=l 

Combining Eqs. (19.11) and (19.12) we get : 


-R 

poo * 11 OO * * r^| 


n w /2 


= 2 • Au; • J3 ™to(0 * 3/(0 


(19.13) 


i=l 


After rearranging terms in Eq. (19.13) we get for the weighting factor An; 

-R 


A w = 


2 • Poo • Uoo ■ a s w • r w • J27={ 2 W w(i) ■ y(i ) 


(19.14a) 


where the number n w of line doublets of the wing is even. Ass uming that the number n v 
of line doublets of the wing is odd, we get : 


An; = t : r (19.146) 

2 ■ Poo * ^oo ■ As w * yi~~i *^w(0 ■ ?/(*") 

Finally, non-symmetric Une doublet strength values <r**(y(i)) that satisfy the mea- 
sured lift, pitching moment, and rolling moment may be determined by using Eqs. (19.5b), 
(19.5c), (19.5d), (19.8a), (19.14a) and (19.14b) . 

The angle of attack correction a t - along the 3/4-chord line of the wing may now be 
computed by applying the wall signature method as the line doublet representation of the 
wing is known (for more detail see Chapter 2) . 

ROLLING MOMENT COEFFICIENT CORRECTION : 

It is now possible to compute a rolling moment coefficient correction due to wind 
tunnel wall interference effects as the angle of attack correction along the 3/4-chord line of 
the wing is known. Figure 50 shows the angle of attack correction distribution a;(0.75, y) 
as a function of the spanwise coordinate y of the wing. Using Fig. 50, we get for the 
infinitesimal rolling moment of a wing span increment d y : 

dR = (-1) • y • - j - J — • a t ‘(0.75, y) ■ qoo ■ dS (19.15) 
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The infinitesimal area d S of the wing span element may be written as : 


d S = c(y ) • d y 


(19.16) 


where c(y) is the local chord of the wing. Combining Eqs. (19.15) and (19.16) and integrat- 
ing the corresponding equation in the spanwise direction, we get for the rolling moment 
due to wall interference effects : 


A R 



d ciXy ) 
d a 


ai(0.75, y) • c(y) • d y 


(19.17) 


The rolling moment coefficient correction Ac? has to be added to the uncorrected 
rolling moment coefficient (see Eq. (19.1)). Therefore we get : 


Ac, = (-1) 


AR 

?oo • S ■ 6/2 


(19.18) 


Finally, combining Eqs. (19.17), (19.18), the rolling moment coefficient correction is ob- 
tained as : 


2 f b / 2 d cr(v) 

Ac, = -= — ■ y ■ — - • a t -(0.75, y) ■ c(y) • dy (19.19a) 

S • b J-b/2 da 

In general, Simpson’s Rule may be applied to Eq. (19.19a) if Ac, has to be determined 
for a given wing. However, it is also possible to simply discretize Eq. (19.19a) as long as 
the local angle of attack correction a, is computed at a large number of equally spaced 
points along the 3/4-chord line of the wing. Assuming that a total number of M equally 
spaced reference points were selected along the complete span of the wing (i.e. b = M ■ Ay 
and dy « Ay = 6/M), we get the approximation : 


M 


Ac, 


-L ■ — ■ y y (i) 

s M 4-* K ’ 

3 - 1 


where 


y(j) = 


d cl(J) 
d a 


ai (0.75, j) • c(j) 


-b 


b 

+ M 


3 2 


(19.196) 


(19.19c) 


The application of Eq. (19.19b) requires an estimate of the local lift curve slope 
dcL(j)/da at the spanwise station of a reference point. This estimate may be obtained by 
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using the Kutta/Joukowski formula. For an infinitesimal wing span increment dy we get 
for the lift dL using the non-symmetric lift distribution : 

d L = poo ■ ■ T**(y) ■ d y (19.20a) 


The local lift coefficient is defined as : 


cl(v) 


dL 

Poo/2 • . c (y) . d y 


Combining Eqs. (19.20a), (19.20b) we then get : 


(19.206) 


c Uy) 


2 • r**(y) 

«oo • c(y) 


(19.21) 


Combining Eq. (19.5a) with Eq. (19.21) and differentiating both sides of the resulting 
equation with respect to angle of attack a we get : 


d cl(p) 


w w (y) 


dr„ 


(19.22) 


da «oo • c(y) d a 

Introducing the abbreviation y(j) = j, we get at the discrete spanwise position j : 

d cl(j) _ 2 • w^(j) d r w 

da Uqo • c(j) d a 

Assuming that the-total lift of the model is charged to the wing, we get from Eq. (19.6c) : 


(19.23) 


L = poo • Uoo • • 6 = c L ■ ■ ula ■ S (19.24) 

Zt 

Differentiating both sides of Eq. (19.21) with respect to the angle of attack a and 
rearranging terms we get : 


d r„ 

d a 


Uoo • S d cl 
2 • 6 ‘ "dlT 


(19.25) 


Thus, combining Eqs. (19.23), (19.25) we get the following approximation of the local 


lift curve slope : 


d c x (j) ^ Ww(j) • 5 d cl 

d a c(j) -6 da 


(19.26) 
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Finally, combining Eqs. (19.19b), (19.26) we get for the rolling moment coefficient 
correction due to wall interference effects : 


A ci » — 


dcL 

da 


M 


M 


^2 * W w(i) * «i(0.75,j) 

i=i 


(19.27) 


where dci,jda is the lift curve slope of the wing in [1/deg], b is the wing span, M is the 
number of reference points in the spanwise direction of the wing, y(j) is the spanwise 
coordinate of reference point wZt(j) is the weighting factor of the non-symmetric 
lift distribution of the wing at reference point “j”, and a,- (0.75, j) is the angle of attack 
correction at reference point “j” in [deg] . An interesting interpretation of Eq. (19.27) 
is possible if we introduce the absolute value of y(j ) and assume that y(j) ^ 0. Then 
Eq. (19.27) becomes : 


Acj 


M 

E 


y(j) |y(i)l 


M |y(j)| 6/2 


/ X dc L 

W w (j) 


da 


oti(0.75,j) 


(19.28) 


Equation (19.28) may be interpreted as the arithmetic mean value of weighted rolling 
moment coefficient contributions along the 3/4-chord of the wing where 

y(j ) 

= sign of the local rolling moment correction 


I y(i) I 


I y(j ) I 

= dimensionless rolling moment arm 


6/2 


{j ) 


dcL 

da 


a,(0.75,j) = local lift coefficient 


Equation (19.27) is used by WICS to determine a rolling moment coefficient correction 
due to wind tu nn el wall interference effects. 
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APPENDIX 20 


SCALE FACTOR LAW 


The minimum of the blockage factor of a wind tunnel model tested in a clean configu- 
ration (i.e. flaps, slats not extended and no landing gear) is more or less identical with the 
solid volume blockage contribution. This approximation is valid if such a model is tested 
in the vicinity of 0.0° angle of attack where solid volume blockage effects axe significantly 
larger that corresponding separation wake blockage effects. 

Solid volume blockage effects are a function of the volume displacement of the wind 
tunnel model. The solid volume blockage factor is proportional to the volume of the test 
article. Therefore, it is possible to develop a scale factor law, that relates the length scale 
of two models of different size but identical geometry to the observed minimum of the solid 
volume blockage. The scale factor law may be written as : 


»(^a) 


Ml) 


Jp2_ 


( 20 . 1 ) 


where e m in is the minimum of the blockage factor as a function of the length scale of the 
model, ip\ is the length scale of the first model, and ip2 is the length scale of the second 
model. 

The following example shows how this scale factor law may be applied to a test of 
two different sized Boeing 7J7 semispan models (see Chapter 3.3 for more detail). Both 
semispan models had identical geometry. The first model was a 8 % scale model, the second 
model was a 14 % scale model. Both models were tested in a cruise configuration. From 
Fig. 23a we estimate that the minimum of the blockage factor computed by WICS (wall 
signature method) and the two-variable method for the 8 % scale model is approximately 
0.0030 . Now it is possible to predict the solid volume blockage factor of the 14 % model 
by applying the scale factor law. Using the blockage factor minimum of the 8 % model 
and the cube of the scale factor ratio we get the prediction : 

14% 


,(14%) 


,( 8 %) 


8% 


= 0.0030 • 5.36 = 0.0160 (20.2) 


Comparing the predicted solid volume blockage factor, i.e. 0.0160, with the estimate 
of the solid volume blockage from Fig. 23b, i.e. 0.0155, we see that both values agree 
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fairly well. Therefore, it can be concluded that the wall pressure measurements used by 
WICS to determine blockage corrections and the corresponding description of the solid 
volume blockage effect by a source/sink pair is sufficiently accurate. The scale factor 
law also provides an additional absolute magnitude check of blockage corrections that are 
computed by using wall pressure measurements. 
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FIGURES 
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REAL-TIME 

LIFT, PITCHING MOMENT, THRUST, 
WALL PRESSURE MEASUREMENT 
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Fig. 1 Basic Elements of Wall Signature Method. 



WIND TUNNEL WALL INTERFERENCE 
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Fig. 2a Definition of Wind Tunnel Wall Interference Effects. 



EMPTY TUNNEL CALIBRATION 



SUPPORT SYSTEM CALIBRATION 



WIND TUNNEL FLOW FIELD 



FREE -AIR FLOW FIELD 

Fig. 2b Singularity Representation of Flow Field. 
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Fig. 2d Wall Interference and Direct Influence of Support System. 
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SINGULARITY REPRESENTATION 
^ OF TEST ARTICLE 
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DYNAMIC PRESSURE CORRECTION 
MACH NUMBER CORRECTION 
ANGLE OF ATTACK CORRECTION 
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Fig. 5b Least Squares Fit of Wall Signature of Fullspan Model. 





SINGULARfTY REPRESENTATION 
OF SUPPORT SYSTEM 
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SUPPORT SYSTEM 


X 
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Fig. 6b Least Squares Fit of Wall Signature of Support System. 



IMAGE PLANE CONFIGURATION 
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Fig. 7 Bipolar Wind Tunnel Approximation of image Plane Configuration. 



IMAGE PLANE CALIBRATION 



Fig. 8 Singularity Representation of Semispan Model. 
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DYNAMIC PRESSURE CORRECTION 
MACH NUMBER CORRECTION 
ANGLE OF ATTACK CORRECTION 
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Fig. 9b Least Squares Fit of Wall Signature of Semispan Model 





Point Source at x = 1 18.0 [ft] , y = 3.0 [ft] , z = 2.0 [ft] 



1050 x [ft] 1340 


Fig. 10 Wall Signature of Point Source as Function of Mach Number. 



Fig. 11a Comparison Bessel Function Solution / Panel Code Solution (WICS). 
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Fig. 11b Comparison Bessel Function Solution / Panel Code Solution (WICS). 
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Fig. 12a Panel Model of Wind Tunnel Geometry / Fullspan Model Configuration. 



Fig. 12b Panel Model of Wind Tunnel Geometry / Semispan Model Configuration. 
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ig. 13b Real-Time Wall Interference Correction System / Semispan Model Configuration. 
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Fig. 14b Pressure / Temperature Measurement during 



WALL PRESSURE PORT 
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Fig. 14 c Application of Matching Conditions to Velocity, Force, and Moment Measurements. 




ROW 8 



LOOKING UPSTREAM 


Fig. 15 Testsection, Semispan Model, and Image Plane Geometry. 
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TS1 15.82 


7J7 8% SEMI-SPAN SINGULARITY 
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*[/<] 
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■■■■ 
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■1 
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-1.000 


SOURCE 
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-2.930 


4 

SOURCE 

■m 

o.coo 

-1.430 


5 

LINE DOUBLET 

119.89 


-3.587 

0.997 

6 

LINE DOUBLET 

120.25 


-2.901 

0.977 

7 

LINE DOUBLET 

120.60 

0.000 

-2.216 

0.934 

8 

LINE DOUBLET 



-1.530 

0.866 

9 

LINE DOUBLET 

121.32 


-0.844 

0.766 

10 

LINE DOUBLET 

121.68 

0.000 

-0.159 

0.619 

11 

LINE DOUBLET 

122.03 

0.000 

0.527 

0.371 


Fig. 16a Singularity Representation of Boeing 7J7 8% Scale Semispan Model. 
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7J7 14% SEMI-SPAN SINGULARITY 
AND REFERENCE LOCATIONS 
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0.934 
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LINE DOUBLET 
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0.000 
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0.866 
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LINE DOUBLET 
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10 

LINE DOUBLET 
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0.000 

2.670 

0.619 

11 

LINE DOUBLET 
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0.000 

3.870 
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Fig. 16b Singularity Representation of Boeing 7J7 14% Scale Semispan Model. 
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1 .0965 





Fig. 17a Wail Pressure Port Location in 12ft Pressure Wind Tunnel. 


ROW 8 ROW 1 



ROW 5 ROW 4 

Test Section Cross Section at x = 120.71 [It] 

(Looking Downstream) 

Fig. 17b Wall Pressure Port Row Definition in 12ft Pressure Wind Tunnel. 
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Wall Pressure Port Number 
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90 


Fig. 19c Wall Signature at Row 3 / 8% Scale Model at 19.82[deg] Angle of Attack. 



Wall Pressure Port Number 


Fig. 19d Wall Signature at Row 6 / 8% Scale Model at 19.82[deg] Angle of Attack. 
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(Row2) 


Wall Pressure Port Number 


Fig. 20a Wall Signature at Row 1 / 14% Scale Model at 9.98[deg] Angle of Attack. 



Woll Pressure Port Number 


Fig. 20b Wall Signature at Row 2 / 14% Scale Model at 9.98[deg] Angle of Attack. 

























Angle of Attock Correction [deg] -|| Angle of Attock Correction [deg] 



-10 0 10 
Uncorrected Angle of Attack [deg] 


ig. 22a Comparison of Angle of Attack Correction / 8% Scale Model. 



Uncorrected Angle of Attock [deg] 

Fig. 22b Comparison of Angle of Attack Correction / 14% Scale Model. 
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Blockoge Foctor [ ] Blockage Factor 



Fig. 23a Comparison of Blockage Factor / 8% Scale Model. 
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Fig. 23b Comparison of Blockage Factor/ 14% Scale Model. 
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DYNAMIC PRESSURE CORRECTION fPSFJ 
Boeing 7J7 Model / 8 Percent Scale 
RUN:SEQ = 0154.0000:041.0000 ; TUNNEL PLANE Y 
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Fig. 24a Singularity Representation and Dynamic Pressure Correction of 8% Scale 
Model at 19.82 [deg] Angle of Attack / Tunnel Plane Y=0.0 . 




DYNAMIC PRESSURE CORRECTION 
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Fig. 24b Singularity Representation and Dynamic Pressure Correction of 14% Scale 
Model at 9.98 [deg] Angle of Attack / Tunnel Plane Y=0.0 . 
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Fig. 24c Singularity Representation and Angle of Attack Correction of 8% Scale 
Model at 19.82 [deg] Angle of Attack / Tunnel Plane Y=0.0 . 




ANGLE OF ATTACK CORRECTION 
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Fig. 24d Singularity Representation and Angle of Attack Correction of 14% Scale 
Model at 9.98 [deg] Angle of Attack /Tunnel Plane Y=0.0 . 




Initial Singularity Representation of Ames Bipod 

No. | Type I x [/<] I y [ft] x [ ft ) Weight 

~l SOURCE 120.5 0.0 -4.0 1.0 

2 SOURCE 120.5 0.0 -3.5 1.0 
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Fig. 25a Singularity Representation of Ames Bipod. 






Fig. 25b Wall Signature of Ames Bipod ; Row 2 . 



Fig. 25c Wall Signature of Ames Bipod ; Row 4 




BLOCKAGE FACTOR [ ] 

AMES-.BIPOD r , 

RUN:SEQ = 0000.0000:011.0000 ; TUNNEL PLANE Y = 0.00 [FT] 
ptot = 4.00 [atm] ; mach = 0.29 ; alpha = 0.00 [deg] 

NO OPTIMIZATION OF SINGULARITY LOCATION 
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Fig. 25d Blockage Factor of Ames Bipod ; Tunnel Plane Y = 0.0 [FT] . 
(No Optimization of Singularity Location) 
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Fig. 25e Blockage Factor of Ames Bipod ; Tunnel Plane Y = 0.0 [FT] . 
(Optimization of Singularity Location) 
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Fig. 25f Blockage Factor of Ames Bipod ; Tunnel Plane X = 120.71 [FT] . 
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Fig. 25g Blockage Factor of Ames Bipod ; Tunnel Plane Z = 0.0 [FT] . 
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Fig. 26 Superposition of Boundary Value Problem Solutions. 
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Fig. 27a Geometry of Rankine Body of Example. 


Pressure Coefficient 


Fig. 27b Panel Representation of Rankine Body Geometry. 
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u(Row4) 


Fig. 29 Panel Model of Rankine Body inside Test Section. 



Fig. 30 Comparison of Wall Signature of Rankine Body with equivalent 
Source-Sink Representation. 
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Fig. 33 Singularity Representation of Rectangular Wing. 



Fig. 34 Comparison of Wall Signature of Wing with 
equivalent Line Doublet Representation. 
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1.33 [ft] I 1.33 [ft] I 1.33 [ft] 





168 


Fig. 35 Comparison of Horseshoe Vortex with Line Doublet. 
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Fig. 36b Shift of Point Doublet Solution. 




Fig. 38 Calculation of Line Doublet Strength using Lift Force and Pitching Moment. 




Fig. 40 Quality Check of Least Squares Fit of Wail Signature. 
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Fig. 41a Linear Interpolation of Wall Pressure Port Calibration (Second Variable) 





High Angle of Attack Sting. 




Fig. 42d Axis of Rotation of Semispan Model. 


Fig. 42e Left Wing Semispan Model. 
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Fig. 43c Pitch a 



Fig. 43e Pitch Rotation of Unit Vector K . 
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Fig. 44 Pitching Moment Lever Arm Calculation. 
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Inclination of Lift and Drag Force. 
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Fig. 45b Inclination of Rolling and Yawing Moment. 
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